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PREFACE 


The  purpose  of  this  report  is  to  give  example  calculations  and  consid¬ 
erations  that  may  be  used  in  the  state-of-the-art  application  of  stone  columns 
for  improving  soils  subject  to  static,  earthquake,  and  post  earthquake-reduced 
strength  loading  conditions  beneath  and  adjacent  to  existing  structures.  This 
report  presents  excess  pore  pressure  considerations  and  the  directly  pertinent 
but  previously  unpublished  appendixes  to  the  following  report  published  by  the 
Federal  Highway  Administration  (FHWA) :  Design  and  Construction  of  Stone  Col¬ 
umns  ,  Vol  I ,  FHWA,  Report  No.  FHWA/RD-83/026 ,  Final  Report,  December  1983, 

210  pp  (by  R.  D.  Barksdale  and  R.  C.  Bachus).  The  appendixes  to  this  report 
were  to  be  published  by  the  FHWA  as  Vol  II,  but  budget  limitations  have  pre¬ 
vented  their  publication;  they  are  given  in  the  present  report  as  Appendixes  A 
through  E. 

The  FHWA  Stone  Column  report  considers  only  briefly  liquefaction  of 
sands  due  to  strong  motion  earthquakes  as  influenced  by  stone  column  construc¬ 
tion.  Therefore,  a  method  is  presented  for  estimating  the  dissipation  of 
excess  pore  pressures  due  to  radial  drainage,  and  practical  construction 
aspects  of  drainage  are  considered. 

This  report  was  prepared  under  Contract  No.  DACW39-85-M-2358  with  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES)  during  the  period  July 
1985  to  December  1986.  The  investigation  was  conducted  under  the  Repair, 
Evaluation,  Maintenance,  and  Rehabilitation  (REMR)  Research  Program  Work 
Unit  32275,  "Remedial  Improvement  for  Liquefiable  Foundations."  Mr.  Arthur  H. 
Walz,  Headquarters,  US  Army  Corps  of  Engineers  (H9DSACE) ,  was  Technica’  Moni¬ 
tor  for  HQUSACE  for  this  work. 

The  REMR  Overview  Committee,  which  approved  this  study,  consisted  ot 
Mr.  John  R.  Mikel  (DAEN-CWO-M) ,  Mr.  Bruce  L.  McCartney  (DAEN-CWH-D) ,  and 
Dr.  Tony  C.  Liu  (DAEN-ECE-D) .  Coordinator  for  the  Directorate  of  Research  and 
Development  was  Mr.  Jesse  A.  Pfeiffer,  Jr.  (DAEN-RDC) ,  and  the  REMR  Program 
Manager  was  Mr.  William  F.  McCleese,  Concrete  Technology  Division,  Structures 
Laboratory,  WES. 

Mr.  Richard  H.  Ledbetter,  Earthquake  Engineering  and  Geophysics  Divi¬ 
sion  (EEGD) ,  Geotechnical  Laboratory  (GL) ,  was  Technical  Monitor  for  the  WES, 
under  the  supervision  of  Dr.  Arley  G.  Franklin,  Chief,  EEGD,  and  under  the 
general  supervision  of  Dr.  William  F.  Marcuson  III,  Chief,  GL.  Appreciation 


is  extended  for  Mr.  Ledbetter’s  help  and  careful  review  of  the  manuscript. 
Special  thanks  also  go  to  Dr.  P.  F.  Hadala,  Assistant  Chief,  GL,  for  his 
thorough  review  of  the  manuscript.  Appreciation  is  expressed  to  the  Federal 
Highway  Administration  for  granting  special  permission  to  allow  the  publica¬ 
tion  of  the  appendixes  to  the  stone  column  report.  Mr.  A.  F.  DiMillio  was 
Project  Manager,  and  Mr.  Jerry  DiMaggio  was  Technical  Monitor  for  the  Federal 
Highway  Administration. 

COL  Dwayne  G.  Lee,  CE,  was  the  Commander  and  Director  during  the  publi¬ 
cation  of  this  report.  Dr.  Robert  W.  Whalin  was  Technical  Director. 


2 


»;  tie'  <**  It-  *. ■»  l.f  t  ,t  «,<  «.» 


t,»  <«>»  »>« 


CONTENTS 


PREFACE. 


CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT . . 


PART  I: 


INTRODUCTION. 


PART  II: 


Stone  Column  Soil  Improvement  Mechanisms — Earthquake 

Loading . . . 

Summary . 

I:  RADIAL  DISSIPATION  OF  EXCESS  PORE  PRESSURES  DEVELOPED 

DURING  EARTHQUAKES . 


Introduction . . . . . 

Liquefaction  Design  Charts  Considering  Radial  Drainage. 
Drain  Permeability . 


PART  III:  SUMMARY  AND  CONCLUSIONS. 


REFERENCES. 


APPENDIX  A:  SELECTED  CONTACTS  FOR  STONE  COLUMNS. 


APPENDIX  B:  LOCAL  BEARING  FAILURE  OF  AN  ISOLATED 
STONE  COLUMN . 


Local  Bearing  Failure . 

Lateral  Bearing  Failure  in  Cohesive  Soil. 

Calculation  of  Limiting  Shear  Force . 

Design  Charts . 


Page 


Design . . . . .  B 1 3 


APPENDIX  C:  EXAMPLE  BEARING  CAPACITY  PROBLEMS. 


Bearing  Capacity  Example  1. 
Bearing  Capacity  Example  2. 


APPENDIX  D:  EXAMPLE  SETTLEMENT  PROBLEMS. 


Settlement  Example  1. 
Settlement  Example  2, 


APPENDIX  E:  EXAMPLE  STABILITY  PROBLEM. 


APPENDIX  F:  EXAMPLE  PROBLEM  OF  PORE  PRESSURE  DISSIPATION  DURING 
AN  EARTHQUAKE . 


APPENDIX  G:  NOTATION. 


ri 

CM 

fell 

K|.<5 

'lV  ’ 

l  I’  •  * 

MV 

;.«.v 

i 

r  * 

Wv 

& 


SfSS 

I  1 

1$ 


m 


CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


horsepower  (550  foot-pounds 
(force)  per  second) 

inches 

kips  (force) 

kips  (force)  per  square  foot 

kips  (mass)  per  cubic  foot 

pounds  (force)  per 
square  foot 

pounds  (force)  per 
square  inch 

pounds  (mass) 

pounds  (mass) 
per  cubic  foot 

square  feet 

tons  (2,000  pounds,  mass) 


745.6999 

2.54 

4.448222 

47.88026 

16018.46 

47.88026 

6.894757 

0.4535924 

16.01846 

0.09290304 

907.1847 


watts 

centimetres 

kllonewtons 

kilopascals 

kilograms  per 
cubic  foot 

pascals 

kilopascals 

kilograms 

kilograms  per 
cubic  metre 

square  metres 

kilograms 
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APPLICATIONS  OF  THE  STATE  OF  THE  ART  OF  STONE  COLUMNS— LIQUEFACTION, 
LOCAL  BEARING  FAILURE,  AND  EXAMPLE  CALCULATIONS 


PART  I :  INTRODUCTION 

1.  Stone  columns  offer  one  method  of  rehabilitating  existing  structures 
to  withstand  both  static  and  earthquake  loadings.  They  can  also  be  placed 
adjacent  to  structures  both  above  and  below  the  water  surface.  Special  tech¬ 
niques  are  available  for  constructing  stone  columns  (or  sand  compaction  piles) 
below  water  (Barksdale  and  Dobson  1983,  Barksdale  1987,  Aboshi  et  al.  1979). 
Stone  columns  can  be  used,  with  careful  consideration  of  seepage  conditions, 
beneath  existing  embankments  and  dams.  Although  they  can  be  installed  through 
existing  dams,  there  Is  a  real  question  of  whether  or  not  they  should  be. 

2.  Construction  of  stone  columns  within  an  earth  dam  can  increase  pore 
pressures  below  the  dam  particularly  near  the  upstream  face  when  the  embank¬ 
ment  is  impervious,  or  when  an  impervious  layer  of  clay  underlies  the  dam.  On 
the  downstream  side  of  the  dam,  flow  is  concentrated  at  the  stone  columns  which 
would  act  as  drains.  Seepage  erosion  could  take  place  at  this  location  of 
finer  particles  into  the  stone  column  drain  under,  for  example,  steady  state 
seepage  conditions.  This  condition  of  seepage  erosion  would  be  further  aggra¬ 
vated  in  regions  where  flow  to  the  stone  column  is  even  further  concentrated, 
such  as  in  a  pervious  layer  near  its  contact  with  an  impervious  layer.  These 
important  aspects  which  have  a  direct  impact  on  dam  safety  certainly  deserve 
careful  consideration  and  further  study. 

3.  For  many  site  and  loading  conditions,  stone  columns  offer  an  excel¬ 
lent  alternative  to  prevent  liquefaction  and  to  improve  stability  and  bearing 
capacity  of  weak  soils.  They  also  significantly  increase  the  rate  at  which 
consolidation  settlement  occurs.  Stone  columns  decrease  the  magnitude  of  set¬ 
tlement,  with  the  maximum  amount  of  improvement  being  about  50  percent  settle¬ 
ment  reduction  for  soft  cohesive  soils  if  about  35  percent  of  the  soil  is 
replaced  by  stone. 

A.  For  both  static  and  earthquake  loadings,  stone  columns  can  be  used 
to  improve  soft  to  firm  clays,  silty  sands,  and  sands  having  silt  and  clay 
layers.  Stone  columns  can  be  constructed  using  either  wet  or  dry  processes. 

i 

I 


Usually,  a  coarse  stone  having  a  variation  In  grain  size  from  0.5-in.*  up  to 
2-  to  3-in.  maximum  size  is  used  in  the  wet  process.  A  comprehensive  descrip¬ 
tion  of  the  design,  construction,  and  general  utilization  of  stone  columns  has 
been  given  by  Barksdale  and  Bachus  (1983)  and  by  Barksdale  and  Dobson  (1983). 

A  general  summary  comparison  including  advantages  and  disadvantages  of  avail¬ 
able  construction  techniques  for  granular  columns  has  been  given  by  Barksdale 
(1987). 

5.  The  purpose  of  this  report  is  to  give  example  calculations  and  con¬ 
siderations  that  may  be  used  in  the  state-of-the-art  application  of  stone 
columns  for  improving  soils  subject  to  static,  earthquake,  and  post¬ 
earthquake-reduced  strength  loading  conditions  beneath  and  adjacent  to  exist¬ 
ing  structures.  This  report  presents  excess  pore  pressure  considerations  and 
the  previously  unpublished  appendixes  to  the  report  entitled  "Design  and  Con¬ 
struction  of  Stone  Columns,  Vol.  I,"  published  by  the  Federal  Highway  Admin¬ 
istration  (FHWA)  (Barksdale  and  Bachus  1983).  These  appendixes  are  directly 
pertinent  to  the  application  of  stone  columns  for  remedial  treatment  of  lique¬ 
fiable  soils. 

6.  Appendixes  A  through  E  of  this  report  presents  the  material  in  the 
same  order  as  referred  to  in  the  FHWA  report.  Appendix  A  gives  addresses  of 
contractors  and  engineers  involved  with  stone  column  construction.  It  is  not 
an  all  inclusive  list  of  such  organizations  and  individuals.  For  the  United 
States,  the  listing  includes  the  contractors  currently  capable  of  constructing 
stone  columns,  and  for  Europe  and  Asia,  the  major  contractors  that  have  the 
capability.  Appendix  B  describes  a  theory  for  estimating  local  bearing  fail¬ 
ure  of  a  single  stone  column  (or  sand  compaction  pile)  subjected  to  a  shear 
loading.  Appendixes  C,  D,  and  E  give  example  problems  illustrating  the  use  of 
stone  columns  for  bearing  capacity,  settlement,  and  stability  problems, 
respectively. 

7.  Appendixes  B  through  E  are  intended  to  be  used  in  conjunction  with 
the  corresponding  theory,  equations,  figures,  and  tables  presented  in  the 
original  FHWA  report  by  Barksdale  and  Bachus  (1983).  Sections  of  the  FHWA 
report  that  present  the  theory  associated  with  Appendixes  B  through  D  of  this 
report  are  given  in  Table  1  as  a  convenient  cross-reference  guide. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 
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Table  1 

Cross-Reference  Between  Appendixes  of  This  Report  and  Theory  Presented 
in  FHWA  Report  "Design  and  Construction  of  Stone  Columns" 


Appendix 


Topic 


Reference  Sections  of 
FHWA  Report  (Barksdale  and  Bachus  1983) 


B 

Local  Bearing  Failure 

Chapter 

VII,  pp. 

153- 

163 

C 

Bearing  Capacity 

Theory: 

Chapter 

III, 

pp.  38-46 

Design: 

Chapter 

VII, 

pp.  142-145 

D 

Settlement  Analysis 

Theory: 

Chapter 

HI, 

pp.  47-76 

Design: 

Chapter 

VII, 

pp.  145-152 

E 

Stability 

Theory: 

Chapter 

III, 

pp.  76-84 

Design: 

Chapter 

VII, 

pp.  152-165 

F 

Pore  Pressure 

Smear: 

Chapter 

III, 

pp.  73-74 

Chapter 

IV, 

pp.  100-105 

Chapter 

VII, 

pp.  150-151 

Construction  Methods 


Chapter  II 


Note:  The  basic  theory  given  in  Chapter  III,  pp.  27-38,  of  the  FHWA  report  is 
applicable  to  Appendixes  B  through  E. 

8.  In  addition,  the  present  report  considers  the  important  concept  of 
the  dissipation  of  earthquake  generated  excess  pore  pressures  due  to  radial 
drainage  into  stone  columns  (or  stone  drains).  The  theoretical  approach 
developed  by  Seed  and  Booker  (1976)  for  evaluating  excess  pore  pressure  dissi¬ 
pation  is  summarized,  and  a  design  example  is  given  in  Appendix  F.  Although 
pore  pressure  dissipation  concepts  are  presented  in  the  present  report. 

Table  l  gives  a  cross-reference  to  appropriate  sections  of  the  FHWA  report 
describing  smear,  which  is  a  very  important  factor  in  radial  drainage. 

Stone  Column  Soil  Improvement  Mechanisms — Earthquake  Loading 

9.  Construction  of  stone  columns  can  under  certain  conditions  improve 
the  ground  when  subjected  to  a  strong  motion  earthquake  due  to  (a)  densifica- 
tion  of  in  situ  granular  soils  surrounding  the  stone  column,  (b)  the  shear 
strength  contributed  by  the  stone  columns,  and  (c)  the  dissipation  of  excess 
pore  pressures  through  radial  drainage  to  the  stone  columns.  The  concepts 
involved  with  these  three  mechanisms  of  improvement  are  briefly  summarized  in 
this  section. 


Densif ication 


10.  The  heavy  vibrator  used  to  install  stone  columns  densifies  clean 
sands  and  to  a  lesser  extent  dirty  sands.  Some  benefit  is  even  derived  from 
installing  stone  columns  in  sandy  silts,  silts,  and  clays  although  the 
improvement  may  take  several  months  to  become  apparent,  and  is  usually  not 
relied  upon  in  design  because  of  the  lack  of  knowledge  concerning  lateral 
stress  behavior  and  strength  gain  with  time. 

11.  If  a  heavy  vibrator  is  used  having  150  hp  or  more,  the  soil  is  sub¬ 
jected  to  a  relatively  severe  condition  of  vibration  during  stone  column  con¬ 
struction.  This  vibration  level  is  probably  equivalent  to  a  moderate  to 
severe  earthquake. 

12.  The  relative  density  of  sands  achieved  when  stone  columns  are  con¬ 
structed  in  the  field  is  dependent  upon  several  factors:  (a)  size  (horse¬ 
power)  and  type  vibrator,  (b)  stone  column  spacing,  and  (c)  characteristics  of 
the  in  situ  soil  including  fines  content,  grain  size,  and  gradation.  The  type 
fines  (i.e.,  amount  of  silt  compared  to  the  amount  of  clay)  also  influences 
performance.  For  usual  stone  column  construction  equipment  and  spacings  of 
about  6  to  9  ft,  relative  densities  of  70  to  85  percent  can  be  usually 
achieved  in  sands  having  fines  contents  less  than  about  5  percent.  Figure  1 
gives  a  generalized  relationship  between  relative  density  and  sand  character¬ 
istics.  This  relationship  can  be  used  as  a  preliminary  guide  in  design. 

Shear  strength  of  stone  column 

13.  Stone  columns  typically  replace  about  20  to  35  percent  of  the  soil 
being  improved.  Also,  the  coarse,  open-graded  stone  used  in  constructing  the 
stone  column  is  quite  dense  and  not  likely  to  undergo  liquefaction.  There¬ 
fore,  the  shear  strength  of  the  stone  column  can  contribute  a  significant 
amount  of  lateral  force  resistance  during  an  earthquake.  Shear  strength  of 
the  stone  column  was,  for  example,  relied  upon  at  the  Santa  Barbara  sewerage 
treatment  plant  (Engelhardt  and  Golding  1975,  Mitchell  and  Huber  1985, 
Barksdale  and  Bachus  1983)  to  carry  the  design  earthquake  loading. 

14.  To  be  able  to  rely  upon  the  stone  column  for  shear  strength  during 
an  earthquake,  some  degree  of  lateral  support  must  be  present  around  the  col¬ 
umn.  Field  shear  tests  are  an  excellent  approach  for  evaluating  the  resis¬ 
tance  of  the  stone  column  and/or  the  composite  soil  mass.  Also,  in  very  soft 
soils  consideration  should  be  given  to  the  possibility  of  local  bearing  fail¬ 
ure  (Appendix  B) . 


Not*:  1  ft  *  •  0.093  »*) 


"5  70 

c 


2  50 1 


%  ; 


Area  Per  Compaction  Probe  (ft  ) 


Figure  1.  Approximate  variation  of  relative 
density  with  tributary  area — for  preliminary 
design  only  (Barksdale  and  Bachus  1983) 


Pore  pressure  dissipation 


13.  If  the  excess  pore  pressures  generated  during  an  earthquake  are 
dissipated  quickly  enough,  liquefaction  will  not  occur.  The  construction  of 
stone  columns  on  relatively  close  centers  provides  vertical  drains  which 
greatly  reduce  the  radial  flow  path  and,  if  properly  constructed,  can  be  quite 
effective  in  dissipating  excess  pore  pressures.  Part  II  of  this  report  con¬ 
siders  in  detail  the  radial  dissipation  of  excess  pore  pressures  generated  dur¬ 
ing  an  earthquake. 


Summary 


16.  In  design  frequently  just  one  of  the  mechanisms  of  stone  column 
improvement  discussed  in  paragraphs  8  through  14  is  considered  to  resist  earth¬ 
quake  effects.  In  the  past,  for  example,  only  the  shear  strength  due  to  the 
stone  column  has  frequently  been  used  in  design.  In  actuality,  two  or  even 
all  three  of  these  mechanisms  will  probably  actively  contribute  to  the  stabil¬ 
ity  of  the  soil  mass  when  subjected  to  seismic  loadings.  The  relative  contri¬ 
bution  of  each  mechanism  depends  upon  many  factors:  (a)  construction  technique 
used,  (b)  specific  site  conditions,  and  (c)  characteristics  of  the  design 
earthquake.  Having  several  different  mechanisms  of  defense  against  failure  is 
in  many  instances  an  important  advantage  over  other  alternative  approaches 
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which  may  rely  on  a  single  mechanism.  Finally,  caution  should  be  exercised  in 
using  stone  columns  having  a  coarse  gradation  in  water  resource  applications 
where  steady  state  or  transient  seepage  conditions  exist. 


PART  II:  RADIAL  DISSIPATION  OF  EXCESS  PORE  PRESSURES 
DEVELOPED  DURING  EARTHQUAKES 

Introduction 

17.  Methods  for  predicting  the  buildup  of  excess  pore  pressures  In  sand 
due  to  earthquakes  have  been  summarized  for  example  by  Seed  (1976)  and  Finn 
(1984).  Seed,  Martin,  and  Lysmer  (1976)  concluded  that  in  homogeneous  coarse 
sands  O^o  >  ^.7  mm),  the  pore  pressures  generated  by  an  earthquake  may  be 
dissipated  because  of  vertical  drainage  as  quickly  as  they  develop.  In  typi¬ 
cal  alluvial  deposits,  however,  numerous  clay  and  silt  layers  which  will  pre¬ 
vent  vertical  drainage  are  likely  to  be  present.  Installation  of  relatively 
closely  spaced  stone  columns  (or  stone  drains)  composed  of  an  open-graded 
stone  results  in  radial  drainage  (Figure  2),  which  is  not  significantly  influ¬ 
enced  by  soil  stratification. 


0 


0  »  Equivalent  dlaaacar  of  area  tributary  to 

a  (too*  coluan 

Figure  2.  Radial  drainage  to  stone  columns 


18.  Seed  and  Booker  (1976)  have  developed  easy  to  use  design  charts  for 
selecting  the  required  spacing  of  stone  columns  (or  stone  drains)  to  dissipate 
earthquake-generated  pore  pressures  in  sand  layers.  The  geometry  associated 
with  radial  drainage  to  stone  columns  is  shown  in  Figure  2.  As  illustrated  in 


the  figure,  important  geometric  variables  are  the  diameter  D*  of  the  con¬ 
structed  stone  column  and  the  equivalent  diameter  D  of  the  area  tributary 

6 

to  the  column.  For  an  equilateral  triangular  pattern  of  stone  columns,  the 
equivalent  diameter  of  the  tributary  area  is 


and  for  a  square  grid  pattern 


D  =  1.05s 
e 


D  =  1.13s 
e 


where  s  is  the  uniform  spacing  of  stone  columns. 

19.  Now  assume  that  the  design  earthquake  can  be  represented  by  a 

series  of  N  uniform  cycles  of  stress  having  a  magnitude  x  over  a 

period  of  time  t^  .  Further,  assume  the  sand  layer  will  liquefy  under  the 

earthquake  loading  xg^  after  being  subjected  for  undrained  conditions  to 

uniform  stress  cycles  of  the  design  earthquake.  Methods  for  evaluating  these 

design  parameters  have  been  described,  for  example,  by  Seed  (1976)  and  Seed 

and  Idriss  (1982);  refer  to  these  references  for  a  complete  discussion  of  how 

N  ,  x  ,  t,  ,  and  N.  can  be  estimated, 
eq  eq  d  l 

20.  During  an  earthquake  the  induced  excess  pore  pressures  vary  with 
both  time  and  position  within  the  zone  tributary  to  a  stone  column.  Now  let 
u  be  the  induced  pore  pressure  in  the  sand  layer  during  the  earthquake  and 
Oq  the  initial  effective  stress  in  the  layer.  The  greatest  pore  pressure 
ratio  r^  =  u/o^  *-n  san^  layer  during  the  earthquake  can  be  calculated 
for  conditions  of  radial  drainage  using  Figures  3  through  6.  The  average 
value  of  pore  pressure  ratio  r^  taken  over  the  area  tributary  to  the  stone 
column  is  given  in  Figures  7  through  10.  To  determine  either  pore  pressure 
ratio,  the  following  dimensionless  time  factor  is  used: 


4k  tJ 

r  =  -r  d— 
ad  n2 

y  i  ,D 
w  v3 


where 


k^  ■  horizontal  permeability  of  the  sand  stratum 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  G) . 
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Figure  3.  Relationship  between  greatest  pore  pressure 
ratio  and  drain  system  parameters  for  N  /N^  *  1 

(adapted  from  Seed  and  and  Booker  1976) 


Figure  5.  Relationship  between  greatest  pore  pressure 
ratio  and  drain  system  parameters  for  N  /N^  =  3 

(adapted  from  Seed  and  Booker  1976) 


Figure  6.  Relationship  between  greatest  pore  pressure 

ratio  and  drain  system  parameters  for  N  /Nr  =  4 

eq  * 

(adapted  from  Seed  and  Booker  1976) 


Ratio  of  Stooa  Coiuan  DUatccr  co  tributary  Araa  DUaaur,  0/0^ 

Figure  8.  Relationship  between  greatest  average  pore 

pressure  ratio  and  drain  system  parameters  for 

N  /N„  “  2  (adapted  from  Seed  and  Booker  1976) 
eq  l 


uniform  stress  cycles  of  the  earthquake 


time  over  which  the 
are  applied 

unit  weight  of  water 


eq 


mvj  ■  coefficient  of  compressibility 
D  -  diameter  of  the  stone  column 

The  terms  in  the  other  ratios  used  in  the  figures  have  been  previously 
defined.  Liquefaction  in  the  sand  layer  would  begin  to  occur  when  r  *  1  , 
using  the  relationships  in  Figures  3  through  6  for  the  greatest  pore  pressure 
ratio.  In  Figures  7  through  10  similar  relationships  are  given  for  the 
average  pore  pressures  developed  in  the  layer.  For  design  a  suitable  safety 
factor  should  be  applied  to  the  pore  pressure  ratio. 

21.  The  theoretical  development  by  Seed  and  Booker  (1976)  of  the  radial 
drainage  theory  just  presented  includes  the  following  assumptions: 
a.  Vertical  drain  is  infinitely  permeable. 

Darcy's  Law  is  valid. 

Flow  of  water  occurs  in  only  the  radial  direction. 

The  radial  coefficient  of  soil  permeability  is  constant. 

The  coefficient  of  volume  compressibility  is  constant. 

The  change  in  bulk  stress  in  the  sand  is  negligible. 


b. 

c. 

d. 

e. 

f. 
&• 


An  Irregular  cyclic  earthquake  loading  can  be  converted  to  an 
equivalent  number  of  uniform  stress  cycles. 

h.  The  Induced  pore  pressure  ratio  r^  as  a  function  of  the 

cyclic  loading  ratio  (N  /N.)  can  be  approximated  using  the 

eq  x 

approach  of  DeAlba,  Chan,  and  Seed  (1975)  for  an  value  of 

0.7  (refer  to  the  DeAlba,  Chan,  and  Seed  paper  for  the  defini¬ 
tion  of  otj  ) . 


Drain  Permeability 


22.  For  the  radial  pore  pressure  dissipation  theory  just  discussed  to 
be  valid,  the  stone  columns  which  act  as  vertical  drains  must  be  sufficiently 
permeable  so  that  a  buildup  in  excess  pore  pressure  does  not  occur  within  the 
column.  Seed  and  Booker  (1976)  concluded  that  a  drain  permeability  200  times 
greater  than  the  adjacent  sand  being  drained  is  sufficient  to  satisfy  this 
requirement.  As  shown  in  Figure  11  for  a  selected  earthquake,  Seed  and 
Booker's  results  indicate  that  important  beneficial  effects  of  stone  columns 
may  be  achieved  even  if  the  permeability  ratio  is  as  small  as  50  to  100.  It 


T  -  ^ 
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Wj  •  k  of  drain 
k  ■  horizontal  k  of  soil 
Earthquake  Characteristics 
Naq/Nt  -  2 
0/0-5 


Figure  11.  Effect  of  permeability  of  stone  column  backfill  on 
rate  of  pore  pressure  dissipation  (adapted  from  Seed  and  Booker 

1976) 

is  unlikely  that  a  stone  gradation  can  be  selected  that  is  both  free  draining 
and  at  the  same  time  satisfies  filter  criteria  developed  to  prevent  erosion  of 
native  soil  and  subsequent  clogging  of  the  drain. 

23.  An  important  limitation  of  relying  on  stone  columns  for  lateral 
drainage  during  an  earthquake  is  that  local  sand  may  be  deposited  during  con¬ 
struction  within  the  voids  of  stone  columns  installed  by  the  conventional 
vibro-replacement  (wet)  technique.  For  example,  at  Santa  Barbara  (Barksdale 
and  Bachus  1983)  where  stone  columns  were  used,  0  to  2  percent  of  the  stone 
column  material  passed  the  0.5-in.  sieve  when  it  was  delivered  to  the  site. 
After  construction,  however,  the  amount  passing  the  0.5-in.  sieve  had  increased 
to  between  11  and  23  percent.  Such  an  increase  in  sand  content  has  also  been 
informally  observed  after  stone  column  construction  at  other  sites. 

24.  Careful  consideration  must  therefore  be  given  during  design  to  the 
possibility  of  an  important  reduction  in  permeability  of  the  stone  column 
during  construction  when  native  sands  are  present.  Further  research  is  needed 
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to  more  clearly  establish  the  variation  of  gradation  with  depth  within  the 
stone  column  and  the  extent  of  reduction  in  stone  column  permeability  due  to 
deposition  of  native  sand. 

Smear  effects 

25.  In  constructing  stone  columns  (or  other  type  vertical  drains)  in 
cohesive  soils ,  a  zone  of  in  situ  soil  becomes  smeared  and  remolded  at  the 
interface  between  the  column  and  adjacent  soil.  Also,  stone  is  forced  into 
the  soil  surrounding  the*  column.  These  factors  all  reduce  the  permeability  of 
a  zone  around  the  outside  of  the  stone  column  and  hence  can  reduce  its  effec¬ 
tiveness  in  dissipating  excess  pore  pressures.  The  combined  effects  of  smear, 
soil  disturbance,  and  intrusion  is  generally  simply  referred  to  as  "smear." 

26.  A  study  by  Cassagrande  and  Poulos  (1968)  has  shown  that  if  a  casing 
is  advanced  without  the  aid  of  jetting,  the  effective  horizontal  permeability 
of  the  soil  can  be  as  much  as  10  times  less  due  to  smear  effects  than  if  jet¬ 
ting  is  used.  Hence,  use  of  a  method  to  construct  stone  columns  (or  sand  com¬ 
paction  piles)  which  uses  jets  to  aid  in  advancing  the  vibroflot  (or  casing) 
should  cause  much  less  smear  especially  where  cohesive  soils  are  present  than 
a  technique  which  does  not  employ  jets. 

Special  construction  techniques 

27.  As  just  discussed,  the  method  of  stone  column  construction  used  may 
have  an  important  effect  on  (a)  the  gradation  of  the  completed  stone  column 
and  (b)  the  amount  of  smear  and  hence  the  effective  radial  permeability.  Two 
possible  techniques  that  can  be  used  to  minimize  one  or  both  of  these  effects 
are  the  "pull-down"  construction  technique  (Barksdale  and  Bachus  1983,  Barks¬ 
dale  and  Dobson  1983)  and  the  sand  compaction  pile  technique  (Barksdale  1987). 

28.  Both  of  these  methods  avoid  the  potential  problem  of  an  important 
change  in  gradation  of  the  stone  column  occurring  during  construction.  If 
side  jets  are  not  used,  however,  smear  along  the  sides  of  the  column  will  be 
increased,  particularly  where  cohesive  soils  are  present. 

29.  Sand  compaction  piles.  The  sand  compaction  pile  method  of  construc¬ 
tion  can  be  used  in  both  sands  and  soft  clays.  This  technique  involves  driving 
a  hollow  casing  (pipe)  down  to  the  desired  elevation  using  a  vertical  vibrator. 
As  the  casing  is  being  driven,  it  is  generally  filled  from  the  top  with  a 
granular  material  using  a  special  hopper  system  to  handle  the  sand.  A  sand  is 
usually  used  in  Japan,  although  a  coarse  stone  can  just  as  readily  be  employed. 
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Indeed,  crushed  stone  and/or  gravel  have  been  used  in  both  Japan  and  Taiwan  in 
this  type  construction. 

30.  The  sand  compaction  pile  method  of  construction  as  presently  prac¬ 
ticed  does  not  use  jets  on  the  sides  of  the  casing.  Driving  the  casing  down 
without  the  aid  of  jets  would  increase  the  amount  of  smear  occurring  in  cohe¬ 
sive  strata.  Smear  from  cohesive  layers  may  even  be  dragged  down  into  adja¬ 
cent  cohesionless  strata.  Therefore,  smear  in  cohesionless  layers  may  become 
important,  particularly  where  these  strata  are  thin.  The  use  of  jets  on  the 
sides  of  the  casing  may  be  possible  in  sand  compaction  pile  construction, 
which  would  minimize  smear  effects. 

31.  Cased  hole  stone  columns.  Special  "pull-down"  type  equipment  can 
also  be  used  to  construct  stone  columns  at  some  sites  (Barksdale  and  Dobson 
1983,  Barksdale  and  Bachus  1983).  In  summary,  the  pull-down  rig  uses  a  stan¬ 
dard  vibroflot  tube  for  constructing  the  stone  columns.  A  smaller  diameter 
feeder  pipe,  however,  is  attached  to  the  vibroflot.  The  granular  stone  column 
backfill  material  which  may  consist  of  either  stone  or  sand  is  fed  to  the  bot¬ 
tom  of  the  vibroflot  through  the  smaller  feeder  tube  as  the  vibroflot  is 
slowly  withdrawn  from  the  ground.  Although  the  rig  is  designed  to  push  the 
vibroflot  and  feeder  tube  down  through  soft  cohesive  soils,  side  jets  and  the 
vibrator  located  at  the  bottom  of  the  vibroflot  can  also  be  used,  for  example, 
in  sands.  Using  this  equipment,  the  hole  is  fully  supported  throughout  con¬ 
struction,  and  sand  can  even  be  used  as  the  backfill  material. 

32.  For  site  conditions  where  the  pull-down  technique  can  be  used,  the 
problem  would  be  eliminated  of  mixing  local  in  situ  sands  with  the  coarse 
backfill  stone  placed  during  stone  column  construction.  A  maximum  stone  size 
of  only  about  1.5  in.  can,  however,  be  used  with  present  pull-down  type  rigs. 
To  minimize  smear  effects,  the  side  jets  on  the  vibroflot  should  be  used 
throughout  construction  of  the  column. 

33.  Summary .  Careful  consideration  should  be  given  during  design  of 
the  potential  reduction  in  the  effectiveness  of  drains  because  of  the  possi¬ 
bility  of  deposition  of  native  sand  in  the  stone  column  and  smear  effects. 

The  magnitude  of  these  adverse  effects  depends  greatly  upon  the  construction 
technique  used.  Extensive  field  evidence  does  show  that  for  static  load 
applications,  stone  columns  act  as  quite  effective  drains.  A  similar 
extensive  experience  record  demonstrating  the  ability  of  stone  columns  to  act 
as  drains  during  strong  motion  earthquakes  has  not  been  developed.  Particular 


caution  should  be  used  where  stone  columns  are  employed  in  and/or  beneath 
water  resource  structures  in  which  important  hydraulic  gradients  exist  over  an 
extended  period  of  time  due  to  either  transient  or  steady  state  seepage  condi¬ 
tions.  Installation  of  stone  columns  can  cause  increased  pressures  and  could 
lead  to  seepage  erosion  into  the  stone  column.  Essentially  no  experience  has 
been  gained  with  the  use  of  stone  columns  under  these  conditions,  and  further 
study  is  certainly  needed. 
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PART  III:  SUMMARY  AND  CONCLUSIONS 
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34.  Design  charts  are  presented  for  estimating  under  earthquake  loading  j 

conditions  the  radial  dissipation  to  stone  columns  (or  stone  drains)  of  excess  ■ 

( 

pore  pressures  generated  in  sand  layers.  For  these  charts  to  be  valid  the  j 

stone  columns  must  have  a  permeability  about  200  times  greater  than  the  sand  | 

being  drained.  The  construction  techniques  used  to  install  stone  columns  can  1 

t 

have  a  signification  effect  on  the  resulting  permeability  of  both  the  stone 
column  and  surrounding  soil. 

35.  During  design  careful  consideration  should  be  given  to  the  possi¬ 
bility  of  a  reduction  in  stone  column  permeability  which  may  occur  using  the 
conventional  "wet"  vibro-replacement  technique.  This  reduction  in  permeabil¬ 
ity  can  be  caused  where  local  sands  are  present  as  a  result  of  their  deposi¬ 
tion  during  construction  into  the  pores  of  the  coarse  stone  backfill.  Possi¬ 
ble  alternatives  to  the  conventional  wet  technique  are  the  pull-down  and  sand 
compaction  pile  techniques  which  both  use  a  cased  hole. 

36.  The  effective  lateral  permeability  of  cohesive  soils  is  reduced 
during  stone  column  construction  because  of  smear.  Smear  occurs  at  the  inter¬ 
face  between  the  stone  column  and  the  surrounding  cohesive  soil.  It  can  sig¬ 
nificantly  reduce  the  effective  lateral  permeability  of  the  cohesive  soil 
caused  by  remolding  and  intrusion  of  soil  into  the  stone.  Field  studies  have 
shown  that  the  effective  lateral  permeability  of  soils  adjacent  to  drains  con¬ 
structed  by  pushing  a  pipe  in  the  ground  without  jetting  are  about  10  times 
less  than  drains  constructed  using  jetting. 

37.  Where  thin  layers  of  sands  and  clays  are  present,  smear  from  the 
cohesive  soils  may  be  dragged  into  the  sand  layers  and  hence  may  also  reduce 
their  effective  permeability.  Further  research  is  needed  to  better  define 
both  smear  effects  and  changes  in  stone  column  gradation  which  occur  during 
construction. 

38.  For  applications  where  stone  columns  are  subjected  to  a  shear  load, 
a  local  bearing  failure  could  occur  involving  punching  of  the  stone  column 
into  the  surrounding  soil.  For  usual  conditions,  local  bearing  failures  are 
most  likely  to  occur  in  cohesive  soils  having  undrained  shear  strengths  less 
than  about  200  to  400  psf.  For  soils  having  lower  shear  strengths  or  when  an 
angle  of  internal  friction  greater  than  about  42  deg  is  used  for  the  stone 

|  column,  the  theory  given  in  Appendix  B  can  be  used  as  a  design  guide.  The 

l 

! 


local  bearing  failure  theory  presented  Is  for  a  single.  Isolated  stone  column. 

39.  Several  examples  are  presented  illustrating  the  use  of  theory  to 
predict  the  behavior  of  ground  that  has  been  improved  using  stone  columns. 
Examples  of  bearing  capacity,  settlement,  stability,  and  dissipation  of  excess 
pore  pressures  due  to  earthquakes  are  presented  in  Appendixes  C,  D,  E,  and  F, 
respectively.  Except  for  the  excess  pore  pressure  dissipation  theory  pre¬ 
sented  in  this  report,  the  theoretical  development  necessary  to  solve  these 
problems  is  given  in  a  previously  published  report  (Barksdale  and  Bachus  1983) . 
Table  1  of  the  present  report  gives  a  cross-reference  between  the  previously 
presented  theory  and  the  illustrative  examples  given  in  the  appendixes. 
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APPENDIX  A:  SELECTED  CONTACTS  FOR  STONE  COLUMNS 
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Pittsburg,  Pennsylvania  15219 
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Kensetsu  Kikai  Chosa  Co.,  Ltd, 
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3- 14-16,  Nishitenma,  Kita-Ku, 
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APPENDIX  B:  LOCAL  BEARING  FAILURE  OF  AN  ISOLATED  STONE  COLUMN 


1.  Stone  columns  are  an  effective  method  for  resisting  rotational  shear 

failures  involving  soft  clays  in  embankments  and  slopes  (DiMaggio  1978f).  For 
a  conventional  slope  stability  analysis,  the  resisting  shear  force  F  devel¬ 
oped  by  the  stone  column  is  determined  by  multiplying  the  effective  normal 
force,  ,  acting  on  the  shear  surface  by  the  tangent  of  the  angle  of  inter¬ 

nal  friction  of  the  stone,  tan  $  .  The  shear  capacity,  F  ,  of  the  stone 

s 

column  can,  under  unfavorable  conditions,  be  limited  by  a  local  bearing  fail- 
ureff  of  the  stone  column  and  cohesive  soil  behind  the  column  as  illustrated 
in  Figures  B1  and  B2 . 

2.  Now  consider  the  behavior  of  an  isolated,  single  stone  column  sur¬ 

rounded  by  a  cohesive  soil.  If  the  shear  force  in  the  stone  column  is  suffi¬ 
ciently  large  compared  to  the  strength  of  the  surrounding  cohesive  soil,  a 
secondary  failure  surface  can  develop  in  the  stone  column  extending  downward 
from  the  circular  arc  failure  surface  (Figure  Bl).  The  resulting  wedge  of 
failed  stone  is  bounded  above  by  the  circular  arc  failure  surface.  The  lower 
failure  surface  develops  within  the  stone  at  an  angle  resulting  in  the  minimum 

resistance  to  sliding  as  defined  by  force  F  .  The  shear  force  F  applied  to 

the  top  causes  the  wedge  (Figure  B2)  to  slide  downward  and  laterally  in  the 
direction  of  movement  of  the  unstable  soil  mass  above.  Sliding  of  the  wedge 
of  stone  is  resisted  by  the  frictional  resistance  of  the  stone  developed  along 
the  bottom  of  the  wedge  and  the  passive  lateral  resistance  of  the  adjacent 
clay.  If  the  passive  resistance  of  the  clay  is  not  sufficient,  the  stone 
wedge  undergoes  a  local  bearing  failure  by  punching  into  the  clay.  If  a  local 
bearing  failure  of  the  clay  occurs  behind  the  stone  column,  the  capacity  of 
the  column  is  limited  by  the  secondary  wedge  failure.  A  local  bearing  failure 

of  the  clay  behind  the  stone  column  has  been  observed  by  Goughnour  (1981)  dur¬ 

ing  a  direct  shear  test  performed  in  the  field  on  a  stone  column.  Reduced 
strength  of  the  composite  mass  was  also  indicated  at  Santa  Barbara  (Engelhardt 
and  Golding  1975)  and  Steel  Bayou  (Ehrgott  1977). 


Note:  Equation,  table,  or  figure  numbers  with  an  asterisk  refer  to  the  fed¬ 

eral  Highway  Administration  Stone  column  report  (Barksdale  and  Bachus 
1983). 

t  References  cited  in  this  Appendix  can  be  found  on  the  last  page, 
ft  Personal  communication,  September  1981,  between  R.  R.  Goughnour  and 
R.  D.  Barksdale. 
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Figure  Bl.  Wedge  type  local  bearing  failure  of  a  stone  column 


Figure  B2.  Local  bearing  capacity  failure  wedge  in 

stone  column 
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Local  Bearing  Failure 

3.  The  limiting  shear  force  that  can  be  applied  if  a  bearing  failure 
controls  can  be  obtained  for  an  isolated  column  by  considering  the  equilibrium 
of  the  wedge  shown  in  Figure  Bl.  This  wedge  together  with  the  forces  acting 
on  it  are  illustrated  in  Figure  B2.  The  following  symbols  used  in  these  fig¬ 
ures,  as  well  as  in  subsequent  derivations,  are  listed  in  the  Notation 
(Appendix  G) : 

W  =  effective  force  of  stone  in  the  wedge 
s 

y  =  effective  (bouyant)  unit  weight  of  stone  in  wedge 
s 

P  =  ultimate  lateral  resistance  of  the  clay  acting  on  the  wedge 
H 

N,T  *  normal  and  shear  force,  respectively,  exerted  on  the  bottom 
surface  of  the  wedge 

W^,F  =  normal  and  shear  force,  respectively,  exerted  on  the  top  surface 
of  the  wedge 

R  =  radius  of  the  stone  column 
D  =  diameter  of  the  stone  column 

4>g  =  angle  of  internal  friction  of  the  stone  column 

a,B  ”  angle  of  inclination  of  the  lower  and  upper  surfaces  of  the 
wedge,  respectively 

A.  The  upper  surface  of  the  wedge  makes  an  angle  B  with  the  horizon- 
talt.  This  upper  surface  coincides  with  the  circular  arc  failure  surface 
(Figure  Bl).  The  lower  surface  of  the  wedge  makes  an  angle  of  a  with  the 
horizontal.  Now  consider  equilibrium  of  the  wedge.  To  develop  the  required 
relationship  for  F  ,  first  sum  forces  acting  on  the  wedge  in  the  vertical 
direction  and  solve  for  the  unknown  normal  force  N  acting  on  the  bottom  of 
the  wedge  obtaining. 


W  +  W_,  cos  B  +  F  sin  B 
s _ N _ 

cos  a  +  tan  d>  sin  a 


where  the  forces  and  angles  are  shown  In  Figure  B2. 


t  R.  R.  Goughnour  of  the  Vibrof lotation  Foundation  Company  has  previously 
developed  a  solution  similar  in  concept  for  the  special  case  of  B  =  0.  His 
solution  handled  lateral  pressure  on  the  column  slightly  differently  than 

this  solution. 
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5.  Now  sum  the  forces  acting  on  the  wedge  in  the  horizontal  direction,  > 

substitute  for  the  unknown  force  N  using  Equation  Bl,  and  solve  for  the 

limiting  force  F  obtaining 

W..  (sin  B  +  A  cos  g)  +  XW  +  P„ 

F  -  - - - r— — — 5 - S  (B2) 

cos  B  -  A  sin  B 

where 

tan  6  cos  a  -  sin  a 

X  _ _ I® _ 

cos  a  +  tan  d>  sin  a 
s 

3- 

W  **  it  (tan  a  -  tan  8)  R  Y 
s  s 

6.  In  the  derivation  of  Equation  B2,  the  effects  of  adjacent  stone 
columns  and  outward,  lateral  spreading  of  the  stone  columns  were  neglected. 

Neglecting  the  effect  of  adjacent  columns  should  introduce  a  factor  of  conser¬ 
vation  in  predicting  the  effect  of  a  local  bearing  failure  (Broms  1964). 

These  effects  are  offset  by  neglecting  lateral  spreading  which  should  be  on 
the  unconservative  side. 


Lateral  Bearing  Failure  in  Cohesive  Soil 


1 


7.  The  ultimate  passive  pressure  developed  by  the  cohesive  soil  as  the 

wedge  pushes  against  it  can  be  calculated  using  the  theory  presented  by  Broms 

(1964)  for  a  single,  laterally  loaded  pile  embedded  in  a  frictionless  soil. 

As  shown  in  Figure  B3,  the  ultimate  lateral  pressure  at  the  surface  is 

taken  to  be  q,  *  2c  with  the  resistance  increasing  linearly  over  a  depth  of 
h 

three  pile  diameters  where  it  reaches  a  maximum  limiting  value  of  q^  =  9c  . 

The  total  depth  beneath  the  surface  h  +  z  (Figure  B4)  is  considered  in 

o 

determining  the  three  pile  diameters.  Near  the  surface,  the  failure  occurs 
because  of  the  upward  flow  of  cohesive  soil  toward  the  surface.  With  increas¬ 
ing  depth  the  failure  becomes  one  of  the  plastic  flow  of  the  soil  from  the 
front  of  the  pile  around  the  sides  (Figure  B3) . 

8.  For  a  single  rough  pile,  plastic  theory  (Poulos  and  Davis  1980; 

Meyerhof  and  Chaplin  1953)  indicates  that  below  a  depth  of  approximately  three 

diameters  the  ultimate  lateral  capacity  is  about  q,  = 

h 


11  tc  12c.  Use  of 


an  ultimate  resistance  of  9c,  however,  is  felt  to  be  prudent  although  it  may 
be  slightly  on  the  conservative  side.  Further,  the  use  of  q^  =  9c  is  rea¬ 
sonable  since  it  is  equal  to  the  end  bearing  capacity  of  deep  piles  embedded 
in  a  cohesive  soil.  The  value  of  q^  =  2c  used  at  the  surface  is  also  real¬ 
istic  since  it  equals  about  40  percent  of  the  bearing  capacity  of  the  clay  in 
the  vertical  direction. 

9.  Now  consider  the  ultimate  lateral  pressure  developed  on  a  wedge  of 
stone  making  an  angle  a  and  6  with  the  horizontal  as  shown  in  Figure  B2. 
Using  the  pressure  distribution  shown  in  Figure  B3,  the  ultimate  passive  pres¬ 
sure  developed  in  the  clay  for  a  depth  (h  +  z^)  S  3D  as  illustrated  in  Fig¬ 
ure  B4  is 

P„  =  Rcfth  +  z  +  R( 1 . 714  +  tan  a)]  (B3) 

H  J  O 

and  for  a  depth  h  +  z  >  3D 

o 

P„  =  36R2  cY  (B4) 

H 


where 

R  =  radius  of  stone  column 
c  »  cohesion  of  soil 
f  =  tan  a  -  tan  8 

h  ”  depth  of  fill  above  the  stone  column 

z  *  depth  of  the  circular  arc  failure  surface  below  the  top  of  the 
stone  column 

The  sign  convention  used  for  a  and  8  is  shown  in  Figure  B4.  Once  a  trial 
circular  arc  failure  surface  has  been  selected,  the  value  of  B  is  known. 

The  angle  a  is  then  determined  to  give  the  minimum  value  of  shear  force  F 
that  can  be  applied  to  the  top  of  the  wedge  before  a  bearing  failure  occurs. 

Calculation  of  Limiting  Shear  Force 


10.  The  limiting  shear  force  F  in  each  column  for  a  given  circular 
arc  sliding  surface  is  calculated  as  follows: 

a.  Determine  the  angle  8  for  a  critical  circle  and  calculate  the 
effective  normal  force,  WN  (Figure  B4)  at  the  point  on  the 
stone  column  where  the  circular  arc  intersects  the  center  of 
the  stone  column  (Figure  Bl). 

Bb 


b.  Select  at  least  three  trial  values  of  the  angle  of  inclination 
a  of  the  lover  surface  of  the  wedge. 

c.  For  each  value  of  a  calculate  the  ultimate  lateral  soil 

resistance,  ,  using  Equation  B3  or  B4  and  a  representative 

value  of  the  undrained  shear  strength  c  of  the  cohesive  soil. 

d.  For  each  value  of  a  ,  calculate  F  for  a  bearing  failure  in 
the  cohesive  soil  using  Equation  B2. 

e.  Plot  the  shear  force  F  obtained  from  Equation  B2  as  a  func¬ 
tion  of  a  and  select  the  minimum  value  of  F  . 

f.  Calculate  the  shear  force  F  that  can  act  on  the  column  if  a 

local  bearing  failure  does  not  develop:  F  =  W  tan  <f> 

N  S 

£.  If  a  local  bearing  failure  of  the  clay  controls,  the  force  cal¬ 
culated  in  Step  e  will  be  less  than  that  calculated  in  Step  f. 
In  the  stability  analysis,  use  the  smaller  of  these  forces  (or 
reduce  the  value  of  $  used  in  design). 

h.  Repeat  the  analysis  for  several  selected  points  along  the  fail¬ 
ure  surface. 


Design  Charts 

11.  Figures  B5  through  Bl5  present  graphically  the  solution  for  local 
bearing  failure  of  a  single,  isolated  stone  column  for  selected  design  parame¬ 
ters.  The  procedure  for  using  the  charts  is  as  follows: 

a.  Select  tentative  design  parameters  and  perform  a  stability 
analysis  for  the  stone  column  improved  ground.  Plot  the  criti¬ 
cal  circle  through  the  stone  columns.  Examine  for  the  possi¬ 
bility  of  local  failure  several  points  along  the  critical 
circle  where  it  intersects  the  center  of  the  columns.  Measure 
the  inclination  B  of  the  circle  (with  the  correct  sign)  and 

the  depth  h  +  z  of  each  point  (Figure  B4) . 
o 

b.  Calculate  the  effective  vertical  force  W  acting  on  the  stone 
column  at  the  depth  under  consideration  by  multiplying  the  ver¬ 
tical  effective  stress  times  the  area  of  the  stone  column. 
First,  calculate  the  effective  body  stress  due  to  the  stone 
column  at  the  selected  point.  Use  the  bouvant  unit  weight  of 
the  stone  below  the  groundwater  table.  Then  calculate  the  ver¬ 
tical  stress  0  due  to  the  embankment  above  the  stone  column 
and  obtain  the  stress  concentration  in  the  column  usine  0 

s 

=  p^o  Equation  8b*.  Add  the  body  stress  to  and  ir.”1tiply 

by  the  area  of  the  stone  column  to  obtain  W 

v 

c.  Using  W  from  Step  b  and  the  design  value(s)  of  <f  and  the 
cohesion  of  the  clay  c  ,  enter  the  appropriate  figure  and 
estimate  the  value  of  the  reduction  factor  n  . 
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Figure  Bll.  Local  bearing  failure  stability  reduction  factors  for  deep 


failure:  -  42°  ,  c  =  400  psf 
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Figure  B12.  Local  bearing  failure  stability  reduction 
factors  for  shallow  failure:  $  =  42°  ,  c  =  200  psf 


Figure  B15.  Local  bearing  failure  stability  reduction  factors  for  j 

deep  failure:  <j>  ■  50°,  c  =  800  psf  I 

s  i 

d.  The  "deep"  charts  should  be  used  when  the  combined  embankment  i 

height  and  stone  column  depth  h  +  z  is  equal  to  or  greater 
than  three  stone  column  diameters;  o?herwise  the  "shallow" 

charts  should  be  used.  1 

12.  The  ratio  q  ,  obtained  from  these  figures,  is  defined  as  follows: 

F  (from  Equation  B2)  _ 

n  =  - -  (B5) 

V?„  tan  d> 

N  Ts 

Physically  n  Is  the  ratio  of  resisting  force  that  is  developed  by  an  iso¬ 
lated  stone  column  if  a  local  bearing  failure  occurs  to  the  force  developed  if 
local  failure  does  not  occur  (i.e.,  the  force  that  conventionally  would  be 
used  in  a  stability  analysis).  Hence,  q  is  the  reduction  factor  indicating 
when  a  local  bearing  failure  may  become  a  problem  for  the  given  geometry  and 
material  properties  used  in  the  design.  Theoretically,  when  q  <  1  ,  local 
bearing  controls  the  maximum  resisting  force  and  moment  that  can  be  developed 
by  the  stone  column.  A  reduction  in  resisting  force  (and  moment)  developed  by 
the  stone  column  would  result  in  a  reduction  if  safety  factor  of  the  slope 
compared  to  that  computed  for  a  general  shear  failure. 


Design 


13.  Full-scale  and  model  direct  shear  tests  indicate  a  local  bearing 


failure  of  at  least  a  single  stone  column  is  possible.  The  analysis  including 
the  design  curves  just  presented  is  for  a  single,  isolated  stone  column.  The 
relatively  close  proximity  of  adjacent  stone  columns  and  lateral  spreading 
greatly  complicate  the  actual  problem  compared  with  an  isolated  column;  cer¬ 
tainly  further  field  and  model  tests  are  needed  in  addition  to  more  refined 
theories.  Nevertheless,  the  design  charts  and  theory  presented  can  be  used  to 
indicate  when  local  bearing  failure  may  be  a  problem.  Further,  the  proposed 
approach  is  useful  as  a  general  guide  in  design  for  selecting  safe  design 
parameters  (4>  ,  n) . 

14.  The  likelihood  of  a  local  bearing  failure  increases  as  the  shear 
strength  of  the  clay  decreases,  and  as  a  greater  angle  of  internal  friction 

and  stress  concentration  factor  n  are  used  in  design.  For  example,  if  a 
♦  of  42°  is  used,  a  local  bearing  might  occur  in  cohesive  soils  having  un¬ 
drained  shear  strengths  less  than  about  400  psf.  (Examine  Figures  B9  through 
B13  for  typical  values  of  8  and  .)  A  local  bearing  failure  could  occur 

in  higher  strength  cohesive  soils  if  ^  values  greater  than  42°  are  used  in 
design.  Therefore,  when  stability  is  being  analyzed  in  very  soft  and  soft 
cohesive  soils,  the  effect  of  a  local  bearing  failure  on  the  overall  slope 
stability  should  be  considered.  Also,  in  firm  and  stiff  soils  such  an  analy¬ 
sis  may  show  that  use  of  higher  values  of  $  may  be  possible  without  under¬ 

going  a  local  bearing  failure. 

15.  Local  bearing  failure  can  be  easily  handled  in  a  slope  stability 

analysis  using  the  concept  of  a  limiting  angle  of  internal  friction  of 

the  stone.  Using  this  simplified  approach,  several  representative  points  are 

selected  along  the  critical  failure  circle(s)  as  determined  by  a  stability 

analysis  on  the  stone  column  improved  ground.  The  effective  vertical  stress 

W  and  inclination  of  the  failure  circle  6  (with  correct  sign)  at  the 
v 

selected  points  are  determined.  Figures  B8  through  B15  can  then  be  used  to 
determine  if  a  local  bearing  failure  might  occur  at  the  selected  points  (and 
the  actual  magnitude  of  the  reduction  in  the  resisting  shear  force  F) .  If  a 
local  failure  is  found  not  to  occur  over  a  significant  portion  of  the  failure 
surface,  the  design  is  satisfactory;  otherwise  consideration  should  be  given 
to  reducing  41  .  Note  that  the  figures  indicate  local  failure  in  general  may 
be  a  problem  only  when  8<0  (i.e.,  near  and  to  the  outside  of  the  toe  of  the 
slope).  Also,  and  perhaps  more  importantly,  the  charts  serve  to  indicate  when 
local  failure  is  not  of  concern.  In  any  case,  past  experience  and  good 


| 

E? 


1 


engineering  judgment  should  be  taken  into  consideration  in  estimating  the 
stability  of  the  slope. 
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APPENDIX  C:  EXAMPLE  BEARING  CAPACITY  PROBLEMS 


Bearing  Capacity  Example  I 

1.  Example  1  illustrates  prediction  of  the  load  due  to  a  wide  fill  that 
can  be  supported  by  stone  column  improved  ground  to  avoid  a  shear  failure  of 
the  stone  columns.  The  specific  problem  is  to  determine  what  height  of  fill 
the  stone  column  improved  ground  can  safely  support.  Both  a  general  shear 
failure  and  a  local  bulging  failure  in  a  deep,  very  soft  clay  layer  (Fig¬ 
ure  Cl)  must  be  considered.  The  subsurface  conditions  and  pertinent  parame¬ 
ters  needed  to  solve  the  problems  are  shown  in  Figure  Cl.  Assume  the  stone 
column  has  an  angle  of  internal  friction  $  of  42  deg,  and  an  equilateral 
triangular  pattern  of  columns  is  used  having  a  spacing  s  =  7  ft  . 


Figure  Cl.  Bearing  capacity  Example  I — Wide  Fill  Over 
Stone  Column  Improved  Clay 

a.  Area  replacement  ratio  a  .  Calculate  ratio  a  from 
_  s  s 

Equation  5b*: 


Note:  Equation,  table,  or  figure  numbers  with  an  asterisk  refer  to  the 

FHWA  Stone  Column  report  (Barksdale  and  Bachus  1983). 


0.907^  =  0.907  (-3y— f)  =  0.227 


tt  ♦  D2  3.14  (3.5  ft)2 


9.62  ft 


A 

A  =  — 
a 

s 


9.62  ft 
0.227 


42.4  ft  (total  area) 


where 


A  »  stone  column  area 
s 

A  =  total  area  within  the  unit  cell 

Stone  column.  Estimate  the  general  ultimate  capacity  of  the 
stone  column  using  Equation  C2  assuming  a  bulging  failure  occurs 
in  the  upper  three  stone  column  diameters  of  depth.  Since  the 
clay  has  a  PI  <  30  and  is  not  classified  as  very  soft 
(c  <  250  psf) ,  use  N  =22  (Chapter  VII,  Barksdale  and  Bachus 
1983t) .  C 


q  .  =  cN  =  0.45  ksf  (22)  =  9.9  ksf 

ult  c 

P  ,  -  q  .  A  =  9.9  ksf  (9.62  ft2)  =  95.2  kips 

ult  ^ult  s 


where 


q^t  =  ultimate  bearing  capacity  of  stone  c 


olumn 


cohesion  of  soil 


=  ultimate  bearing  capacity  factor  of  stone  column 
In  the  above  expressions,  the  stress  in  the  stone  column  at 


ultimate  is 


0  =  q  ,  =  cN 

s  ult  c 


Deep  bulging.  Now  check  for  the  possibility  of  a  bulging  fail¬ 
ure  in  the  very  soft  clay  stratum  located  at  a  depth  of  20  ft. 

As  discussed  in  Appendix  B,  the  ultimate  lateral  stress  which  an 
isolated  stone  column  can  develop  is  approximately  equal  to 
a  9c  =  9(0.2  ksf)  =K1.8  ksf  since  the  weak  stratum  is  greater 
than  3D  below  the  surface.  From  Equation  9*  the  ultimate  stress 
the  stone  column  can  carry  is  then 


(1  +  sin  ) 
~  _ s 

^ult  °3  (1  —  sin  $  ) 

s 

1  -  =  9.07  ksf 


=  1.8  ksf  (5.04) 


t  References  cited  in  this  Appendix  can  be  found  on  the  last  page. 
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where  equals  minor  principle  stress.  Since  the  ultimate 

stress  the  stone  column  can  carry  considering  a  deep  bulging 
failure  in  the  very  soft  layer  is  slightly  less  than  for  a  fail¬ 
ure  at  the  surface,  the  very  soft  deep  stratum  controls. 

Cohesive  soil.  The  maximum  ultimate  stress  that  the  clay  sur¬ 
rounding  the  stone  column  can  take  is  0  =  5c  =  5(0.450  ksf) 

c 

■  2.24  ksf.  However,  the  total  load  applied  to  the  unit  cell 
must  also  not  overload  the  clay.  Assuming  a  stress  concentra¬ 
tion  factor  n  =  3  ,  from  Equations  8a*  and  8b* 


s  [1  +  (n  -  1)  a  ]  [1  +  (3  -  1)  0.227] 

s 


2.06 


V 

Then 


c  fl  +  (n  -  1)  a  ]  [1  +  (3  -  1)  0.227] 

s 


0.688 


0  S  p  0  =  p 
c  c  c 


0c  <  uc0  =  0.688  f^orf)  =  3.0  ksf 


(8b)* 

(8a)* 

(C3) 


where 

*  average  stress  acting  over  soil  in  unit  cell 

*  ratio  of  stress  in  cohesive  soil  to  average  stress 

og  =  average  stress  acting  over  stone  column 

p  =  ratio  of  stress  in  stone  column  to  average  stress 
s 

Since  3.0  ksf  is  greater  than  5c  =  2.25  ksf,  =  5c 

■  2.25  ksf  is  the  ultimate  stress  the  clay  can  carry. 

Allowable  fill  loading.  The  ultimate  loading  that  can  be 
applied  over  the  unit  cell  area  well  within  the  fill  area  is 

P  ,  =  0  A  +  0  A  =  (9.07  ksf) (9.62  ft2) 

ult  s  s  c  c 

+  (2.25  ksf) (32 .8  ft2)  (C4) 

P  ,  =  161  kips 

ult 

Using  a  safety  factor  of  2.0  the  allowable  loading  is  P 

<x  x.  X 

*  161  kips/2  =  80.5  kips.  The  height  of  embankment  H'  ,  that 
will  apply  the  safe  loading  to  the  unit  cell  is  H'  = 


a 


Hence 


H' 

H' 


p  /  (y^^A) 

all/Uwet  ’ 

15.1  ft 


80,5  kips 


(0.125  kef  x  42.4  ft  ) 


(C5) 


Commentary.  Settlement  of  the  fill  would  be  significant  and 
should  be  calculated.  Also,  the  stability  at  the  edge  of  the 
fill  should  be  checked  using  a  circular  arc  analysis.  In  this 
example  the  very  soft  clay  layer  at  a  depth  of  20  ft  controls 
the  load  that  can  be  applied  to  the  stone  columns.  Use  of  an 
ultimate  lateral  stress  of  9c  acting  on  the  stone  columns  should 
give  a  conservative,  but  realistic,  estimate  of  the  ultimate 
resistance  to  bulging  that  can  be  developed  (refer  to  Appendix  B 
for  a  more  indepth  consideration  of  this  aspect). 

Using  =  9c  as  the  limiting  lateral  pressure  that  the  soil 

can  withstand,  the  ultimate  load  that  a  stone  column  can  carry 
would  for  4>g  =  42°  be  equal  for  depths  greater  than  3D  to 

q"  , .  =  9c  (1  +  sin  $  ) / ( 1  -  sin  4>  )  =  9c  (5.04)  =  45c  or  N 
ult  s  s  c 

=  45  ,  which  indicates  that  a  limiting  value  of  N  exists  at  a 

deep  depth. 

Since  the  fill  is  wide,  the  stress  on  the  stone  column  does  not 
decrease  with  depth  because  of  lateral  spreading  of  stress.  If 
a  narrow  group  of  stone  columns  had  been  used,  the  stress  would, 
however,  decrease  with  depth;  this  could  be  taken  into  account 
to  determine  the  increased  stress  that  could  be  applied  at  the 
surface  compared  with  the  level  of  the  very  soft  clay  stratum 
which  controlled. 

Finally,  Vesic  (1972)  cavity  expansion  theory  could  also  have 
been  used  to  determine  the  ultimate  capacity  of  the  stone  column 
in  the  weak  stratum.  Since  the  clay  is  very  soft  and  has  a 
PI  >  30,  E  *  5c  is  used  to  calculate  a  Rigidity  Index,  I 
(Equation  13*)  of  1.72  for  v  =  0.45  .  In  this  analysis  let  q 
equal  the  total  lateral  stress  acting  at  the  center  of  the  soft 
layer.  Nonlinear  finite  element  analyses  indicate  the  lateral 
pressure  due  to  the  applied  surface  loading  o  can  be  conser¬ 
vatively  approximated  as  0,4o  : 


q  =  K  yz  +  0.4o  =  0.75(24  ft)(0.1  kef)  +  0.4  (2.25  ksf)  (C6) 

o  c 


q  -  2.7  ksf 


Now  F'  *  Jin  I  +  1  =  £n  1.72  +  1  =  1.54  for  p  =  0  ,  Fig- 
c  r  Yc  6 

ure  C2  (16*),  and  no  volume  change.  Then  the  ultimate  load  the 

stone  column  can  carry  is 


i  -  — *1  r~ 2  5  •*  o'* 

Mrs.  Clean  Sand 

Figure  C3.  Bearing  capacity  Example  2 — Square  Group 


3.  From  Figure  C4  (27*)  the  stress  concentration  factor  in  the  stone 
column  improved  ground  is  about  2.0.  The  bearing  capacity  calculations  are  as 
follows : 

a.  Calculate  the  area  replacement  ratio,  a 
-  s 

A  =  ^  (2.5  ft)2  *  10  =  49.1  ft2  (C7) 

s  4 


A  =  13.5  ft  x  10.5  ft 


a 

s 


49.1  ft2 
141.8  ft2 


=  141.8 

0.346 


ft' 


(C8) 

(C9) 


b.  Determine  the  stress  concentration  ratio  in  the  stone  column 
from  Equation  8b*  or  Figure  C5  (68*) 


Strata  Concaotratlon  Factor 


s:| 


St  res*  in  Clay  and  Stona 


«,*#,«  •  *c« 

mc«  1/11+  (n-i)  a^I 


Strcaa  Concentration  Factor,  n 

Figure  C5.  Variation  of  stress  in  clay  with  stress  concentia- 
tion  factor  and  area  ratio  (Barksdale  and  Bachus  1983) 


s  "  [1  +  (n  -  1)  a  ]  ~  [1  +  (2  -  1)(0.346)] 
s 


1.49  (8b 


Calculate  the  composite  shear  strength  within  the  stone  column 
group  (Equations  16a*  and  J6b*)  and  related  parameters  as 
follows ; 

[tan  d> ]  =  u  tan  *  (a  )  =  (1.49)  (tan  42°)  (16a 

avg  s  s  s 

x  (0.346)  -  0.464 

1  ^ 

<f>avg  -  tan"  (0.464)  -  24.9°  and  8  =  45  +  (Cl 

2 

tan  8  =  1.566  and  tan  6  “  2.454  (Cl 

c  -  c  •  (1  -  a  )  =  1  ksf  (1  -  0.346)  *=  0.654  ksf  (16b 


.V.V.V-V.*4 


Using  Vesic  cavity  expansion  theory,  calculate  the  ultimate  lat¬ 
eral  stress  a ^  in  the  clay  surrounding  the  stone  column  group. 
Since  the  clay  is  stiff,  has  no  organics,  and  has  a  PI  =  30,  use 
E  =  1 lc  for  calculating  the  Rigidity  Index,  I  .  The  average 

diameter  of  the  foundation  is  B  =  V4A/3. 14  =  13.4  ft.  The 
depth  of  the  failure  wedge  is  then  (Figure  C3)  B  tan  B  +  3  ft 
=  (13.4  ft) (1.566)  +  3  ft  =  24  ft.  The  initial  lateral  stress 
in  the  stiff  silty  clay  surrounding  the  stone  columns  will  be 
used  as  a  conservative  estimate  of  the  mean  stress  q  (Equa¬ 
tion  12*),  for  use  in  the  cavity  expansion  theory.  The  stiff 
silty  clay  is  known  to  be  normally  consolidated .  Therefore, 
from  Lambe  and  Whitman  (1969),  =  0.6  for  the  surrounding 

silty  clay,  and  q  s  0.6  (13.5  ft  *  0.115  kef)  =  0.931  ksf.  Now 
calculate  the  Rigidity  Index  (Equation  13*). 


I 

r 


_  E _  =  _ lk _ 

2(1  +  v)  (c  +  q  tan  <Ji)  2(1  +  0.45)  (c  7  q  tan  0°) 


(13)* 


giving  1^  =  3.79  .  From  F^  =  £n  I  +  1  for  <p  =  0  and  Fig¬ 
ure  C2  (16*),  F'  =  2.33  and  F'  =  1.0  .  Then  calculate  the 
c  q 

ultimate  lateral  stress  which  can  be  developed  by  the  surround¬ 
ing  silty  clay: 

a,  -  cF'  +  qF'  -  1  ksf  (2.33)  +  0.931  ksf  (1.0)  (12)* 

3  c  q 

-  3.26  ksf 

Calculate  the  ultimate  vertical  stress  and  load  that  can  be 
applied  over  the  rigid  foundation  (refer  to  Equation  19*): 

q  ,  =  o,  =  ct~  tan^  6  +  2c  tan  B  =  3.26  ksf  (2.454) 

ult  1  3  avg 

+  2(0.654  ksf) ( 1 . 566)  (19)* 

^ult  **  8.0  ksf  +  2.0  ksf  =  10.0  ksf 

The  ultimate  load  that  can  be  carried  by  the  foundation  Is  P  .  _ 

2  ult 

■  q  ^  *  A  «  10.0  ksf  (141.8  ft  )  «=  1,418  kips.  Using  a  safety 

factor  of  2.0,  the  foundation  can  carry  P  =  1^18  kips/2.0 

*  709  kips.  This  amounts  to  70.9  kips  (or  35.5  tons)  per  stone 
column  if  the  silty  clay  is  assumed  not  to  carry  any  of  the 
load.  This  level  of  loading  is  reasonable  for  a  foundation 
where  settlement  is  of  concern.  Table  Cl  (12*). 

Settlement  of  course  would  control  the  design.  A  total  load  on 
the  group  of  709  kips  would  be  used  for  a  first  settlement  esti¬ 
mate.  For  this  loading  the  average  stress  applied  to  the  foun- 

2 

dation  is  o  =  709  kips/141.8  ft  =5  ksf.  The  probable 


» 


distribution  of  stress  between  the  stone  and  soil  for  n  =  2 
would  be 

a  =  u  o  =  0.743  (5  ksf)  =  3.7  ksf  (C12) 

c  c 

a  =  u  a  =  1.49  (5  ksf)  *=  7.45  ksf  (C13) 

s  s 

Since  the  ultimate  stress  of  the  stiff  clay  is  about  6.2c 
=  6.2  ksf,  the  stress  level  in  the  clay  is  not  excessive.  Using 
the  proposed  design,  the  ratio  of  the  settlement  of  the  treated 
to  unimproved  ground  would  be  approximately  St/S  £  V*c  =  0.74 

(refer  to  Equations  20*,  21*,  and  22*  and  Figure  C6  (19*)). 

Thus  for  the  conditions  analyzed,  reduction  in  settlement  on  the 
order  of  25  percent  would  be  expected.  For  the  given  site  con¬ 
ditions,  use  of  a  larger  footing  (without  stone  columns)  should 
also  be  evaluated  considering  magnitude  of  settlements  and  the 
economics  of  the  designs. 

In  general  for  the  wet  method,  a  stone  column  spacing  less  than 
5  ft  is  not  recommended;  Example  Problem  2  would  therefore  be  an 
exception  because  of  the  presence  of  the  stiff,  silty  clay. 


Table  Cl 

Approximate  Range  In  Design  Loads  Used  In  Practice  for 
Stone  Columns  (Barksdale  and  Bachus  1983) 


Soil  Type 

Approximate  Design 

Load ,  tons 

Foundation  Design ^ 

Stability^ 

(2) 

1.  Cohesive  soilv 

15-30 

20-45 

400  psf  <  c  <  600  psf 

600  psf  <  c  <  1,000  psf 

25-45 

30-60 

c  >  1,000  psf 

35-60 

o 

r— 

1 

O 

2.  Cohesionless  Soil 

20-180 

— 

(see  Note  1) 

Notes : 


1.  In  general,  when  stone  column  loads  are  given,  all  the  applied 
load  Is  considered  carried  by  the  stone  column. 

2.  Typical  design  loads  for  foundations  on  cohesive  soils  are  15  to 
30  tons . 

3.  Unit  Conversions:  1  psf  *  47.9  N/m 
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APPENDIX  D:  EXAMPLE  SETTLEMENT  PROBLEMS 

Settlement  Example  1 

1.  Settlement  Example  1  illustrates  calculating  settlements  of  a  soft 
clay  reinforced  with  stone  columns  and  loaded  by  a  wide  fill.  The  calculation 
of  the  load  carrying  capacity  of  stone  column  improved  ground  for  a  problem 
similar  to  this  was  illustrated  by  Example  1  in  Appendix  C.  In  the  present 
example,  primary  consolidation  settlements  are  calculated  using  both  the  equi¬ 
librium  and  finite  element  methods.  Secondary  settlements  are  also  calcu¬ 
lated.  The  problem  is  illustrated  in  Figure  Dl.  The  site  consists  of  20  ft 
of  gray,  soft  silty  clay  overlying  a  firm  to  dense  sand.  The  groundwater 
table  is  at  the  surface.  An  equilateral  triangular  pattern  of  stone  columns 


/ 

/  Embankment  (y  •  120  pcf ) 
f  y  wet  r  J 

/  /  /  /  A 


*  Sand  Blanket  (y  “  108  pcf) 

wet  •  • 


Figure  Dl.  Settlement  Example  1 — Wide  Fill  Over 
Stone  Column  Improved  Silty  Clay 


Note:  Equation,  table,  or  figure  numbers  with  an  asterisk  refer  to  the  Fed¬ 

eral  Highway  Administration  stone  column  report  (Barksdale  and  Bachus 
19831-). 

t  References  cited  in  this  Appendix  can  be  found  on  the  last  page. 
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is  used  having  a  spacing  of  6.5  ft.  The  diameter  of  the  stone  column  is 
estimated  from  Table  D1  (13*)  to  be  3.5  ft.  A  2.5-ft  sand  blanket  is  to  be 
placed  over  the  soft  silty  clay  for  a  working  platform  and  drainage  blanket. 
Equilibrium  method 

2.  The  average  stress  a  exerted  by  the  2.5-ft  sand  blanket  and 
12.5-ft  structural  fill  on  the  top  of  the  stone  columns  is  o  =  12.5  ft 
x  120  pcf  +  2.5  ft  x  108  pcf  *  1,770  psf.  The  area  replacement  ratio  ag 
from  Equation  5b*  is  for  an  equilateral,  triangular  stone  column  pattern: 


a  =  0.907 
s 


0.907 


0.263 


(5b)* 


Assume  for  the  settlement  analysis  the  stress  concentration  factor  n  to  be 

✓ 

5.0.  Then  the  stress  concentration  ratio  uc  in  the  clay  is  from  Equa¬ 
tion  8a*  or  Figure  C5  (68*): 


=  _ 1 _ 

[ 1  +  (n  -  1 )  a  ] 
s 


0.487 


(8a)* 


The  initial  effective  stress  oq  at  the  center  of  the  silty  clay  layer  is 

o^  =  10  ft  x  (95  pcf  -  62.4  pcf)  =  326  psf  (Dl) 

The  primary  consolidation  settlement  in  the  clay  layer  from  one-dimensional 
consolidation  theory  is  from  Equation  20*: 


a  +  a 

c  _  _ c  ,  I  o  c 

st  1  +  e  lo8lo' 


(20)’ 


\1  +  2.0/ 


31.4  in. 


log 


10 


326  psf  +  (1770  psf) (0.487) 
326  psf 


(20  ft  x  12  in. /ft) 


where 


S^.  *  primary  consolidation  settlement  occurring  over  a  distance  H  of 
stone  column  treated  ground 
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compression  index  from  one-dimensional  consolidation  test 
initial  void  ratio 


C  - 

c 

e  « 

o 

0q  »  average  initial  effective  stress  in  the  clay  layer 

o  =  change  in  stress  in  the  clay  layer  due  to  the  externally  applied 
loading,  Equation  8a* 

H  *  vertical  height  of  stone  column  treated  ground  over  which 
settlements  are  being  calculated 

The  estimated  primary  consolidation  settlement  of  the  stone  column  improved 
silty  clay  layer  is  thus  31  in.  following  the  equilibrium  method.  For  compar¬ 
ison,  the  settlement  in  the  silty  clay  layer  if  not  improved  with  stone  col¬ 
umns  would  be  S  *  45.2  in. 

3.  Note  how  simple  the  equilibrium  method  is  to  apply  to  a  problem. 

The  "trick,”  of  course,  is  to  estimate  the  correct  value  of  stress  concentra¬ 
tion  factor  n  to  use  in  the  analysis.  In  this  problem,  the  fill  was  wide 
and  no  dissipation  laterally  of  stress  with  depth  occurs.  The  next  settlement 
example  shows  how  both  the  equilibrium  and  the  finite  element  methods  can  be 
applied  to  a  problem  where  the  applied  stress  decreases  with  depth. 

Nonlinear  finite  element  method 

4.  Since  the  clay  is  soft  and  quite  compressible  use  the  nonlinear 
finite  element  method  of  analysis.  First,  calculate  the  modulus  of  elasticity 
E^  of  the  clay  for  the  approximate  stress  range  of  interest.  The  initial 
average  stress  in  the  clay  from  the  equilibrium  method  is  =  326  psf.  The 
change  in  stress  in  the  clay  due  to  the  embankment  loading  is  =  p^o 

=  0.487  (1,770  psf)  =  862  psf.  Using  Table  D2  (9*)  and  experience  as  a  guide, 
the  drained  Poisson's  ratio  of  the  clay  is  assumed  to  be  0.42.  From  Equa¬ 
tion  47*  the  modulus  of  elasticity  of  the  clay  for  the  applicable  stress  range 
is 


E  ** 
c 


X 


(1  +  v) (1  -  2v)  (1  +  e  )a 

o  va 


0.435(1  -  v)  C 

(326  psf  +  862  psf) 
2 


c 

■] 


2,292 


(1  +  0.42) (1  -  2  x  0.42) ( 1  +  2.0) 
0.435(1  -  0.42) (0. 70) 

psf  =  15.9  psi 


(47)* 


where 

E  =  drained  modulus  of  elasticity  (for  a  stress  path  along  the  K 
C  line) 

v  =  Poisson's  ratio  (drained) 


a  =  average  of  initial  and  final  stress  state  applied  in  the  field 
(vertical  stress) 


u1  a,*  ■.  i 


k*  «>  «  <  ».*  i  •  «.»  #,*  9.*  ju*  i.f  h* 


Table  D2 


Typical  Poisson's  Ratio  Values  of  Clay  for 


Drained  Loading  (Poulos  1975) 


Soil  Consistenc 


e  -,(2,3) 


Poisson's  Ratio' 


Very  Soft  to  Soft 


0.35  -  0.45 


Firm  to  Stiff 


0.30  -  0.34 


Stiff  Overconsolidated  Clays 


0.1  -  0.29 


Notes : 


1.  For  undrained  loading  use  0.45. 

2.  For  normally  consolidated  clays. 

3.  For  very  soft  to  soft  clays  a  value  of 
0.40-0.45  is  recommended  for  calculating 
E  for  nonlinear  finite  element  settle- 
ment  analyses  of  stone  column  improved 
ground;  for  firm  to  stiff  use  at  least 

v  -  0.35  . 


Note  that  the  value  of  Poisson's  ratio  selected  has  a  significant  effect  on 


the  calculated  value  of  E  ;  larger  values  of  v  give  smaller  values  of 

c 


5.  The  stone  column  length  to  diameter  ratio  in  the  soft  clay  is  L/D 
*  20  ft/3.5  ft  -  5.7  .  The  average  applied  pressure  o  due  to  the  embankment 
is  a  *  1,770  psf  «  12.3  psi.  Interpolating  from  Figures  32*  and  33*  (a 


0.25)  for  a  soft  boundary  condition  (E,  »  12  psi),  the  settlement  of  a  stone 

D 


column  of  length  L  *  20  ft  and  diameter  D  =  3.5  (L/D  =  5.7)  is  21  in.  In 


interpolating  between  figures  for  different  S^/L  values,  work  in  terms  of 


settlement  when  the  length  varies.  Hence,  calculate  the  actual  settlement 


L  (S^/L)  for  each  value  of  L/D  and  before  interpolating  for 


settlement  using  a  graph.  The  "best"  settlement  estimate  is  the  average  of 
the  finite  element  and  incremental  methods:  S  =  (21  in.  +  31  in.)/2  =  26  in. 


The  estimated  reduction  in  settlement  due  to  stone  column  improvement  is  then 


St/S  -  26.0  in. /45 . 2  -  0.575 


Time  rate  of  settlement 


6.  Determine  the  magnitude  of  primary  consolidation  settlement  after 
2  months  assuming  instantaneous  construction  (Leonards  1962).  (Methods  for 
handling  construction  over  a  finite  time  interval  are  also  presented  by 


Leonards.)  The  silty  clay  has  a  vertical  coefficient  of  consolidation  C  of 


Ss 


.is 


P 


8 


m 


n 


v* 


0.05  ft  /day.  Based  on  a  detailed  study  of  the  strata,  the  horizontal  permea¬ 
bility  is  estimated  to  be  three  times  the  vertical  permeability.  Then  from 
Equation  49* 

Cvr  =  Cv  (ir)  =  (0-°5  ft2/dayH3)  =  °-15  ft2/ day  (49)* 

^  v ' 

where  C  =  horizontal  coefficient  of  consolidation, 
vr 

7.  Assume  the  reduced  drain  diameter  D'  to  account  for  smear  is 
one-fifth  the  constructed  stone  column  diameter.  For  an  equilateral  triangu¬ 
lar  stone  column  spacing  s  of  6.5  ft,  the  equivalent  diameter  of  the 

unit  cell  is 


1.05s  =  1.05  (6.5  ft)  =  6.83  ft 


6.83  ft 
3.5  ft/5 


=  9.76 


(Figure  45)* 


where  n*  equals  ratio  of  unit  cell  radius  r  to  the  radius  of  drain  r 

e  w 

(stone  column  radius  less  smear  zone  thickness).  The  dimensionless  vertical 
and  horizontal  time  factors  are  then 

T  =  — - -  =  (0.05  ft2/day)  (-2  X  3-1-  da.^-)-  =  0.031  (27)* 

Z  (H/N)  (20  ft/2) 


T  .  Ssl  .  (0.15  ft2/day>  *■ -3j- dV2>  .  0.  ,99 


(6.83  ft)' 


where 


T  =  time  factor  for  vertical  direction 
z 

C  =  coefficient  of  consolidation  in  vertical  direction 
v 

t  =  elapsed  time  of  consolidation 
H  =  thickness  of  cohesive  layer 

N  =  number  of  permeable  drainage  surfaces  at  the  top  and/or  bottom  of 
the  layer  (N  =  1  or  2) 


=  time  factor  for  radial  drainage 
Dg  *  equivalent  diameter  of  unit  cell 

From  Figure  D2  (42*)  the  degree  of  consolidation  in  the  vertical  direction  is 
”  0.12  and  from  Figure  D3  (43*)  the  degree  of  consolidation  in  the  radial 
direction  is  =  0.64  .  The  combined  degree  of  consolidation.  Equation  25* 
is 


U  =  1  -  (1  -  U  )(1  -  U  )  =  1  -  (1  -  0. 12) (1  -  0.64)  =  0.68  (25)* 

z  r 

where  U  =  average  degree  of  consolidation  of  the  cohesive  layer  considering 
both  vertical  and  radial  drainage.  An  important  portion  of  the  total  primary 
consolidation  settlement  occurs  before  2  months  and  equals  =  26.0  in. 
(0.68)  *  18  in.  For  this  example  problem  having  stone  columns,  vertical 
drainage  of  the  soil  had  little  effect  on  the  time  rate  of  primary  settlement 
because  of  the  higher  radial  coefficient  of  consolidation  and  smaller  radial 
drainage  path  to  the  vertical  drains.  For  comparison,  if  stone  columns  had 
not  been  used,  primary  consolidation  settlement  would  have  been  only 
12  percent  complete,  with  the  primary  settlement  at  the  end  of  2  months  being 
only  about  3  in. 

Secondary  compression  settlement 

8.  Estimate  the  magnitude  of  secondary  compression  settlement  that 

would  be  expected  to  occur  5  years  after  construction.  Assume  secondary 

compression  begins  at  the  time  for  90  percent  primary  consolidation.  Neglect 

the  effects  of  vertical  drainage  which  were  shown  above  to  be  small.  The 

radial  time  factor  for  90  percent  primary  consolidation  for  n*  =  9.76  is 

T^  »  0.47  from  Figure  D3  (43*).  From  Equation  28*  the  time  for  90  percent 

primary  consolidation  is  t  =  T  (D  )^/C  =  0.47(6.83  ft)^/(0.15  ft^/day) 

r  e  vr 

=  146  days  after  construction.  The  secondary  compression  settlement  is  then 

AS  =  CaH  log1Q  (t2/tx)  (30)* 

AS  *  0.005(240  in.)  log^  [5(365  days)/146  days]  =  1.3  in. 

where 


AS  “  secondary  compression  of  the  layer 


i  ivt*  -r-w  v/w  1 


/^r.  AT.  fvi  *n  WLAXaXaiiaTLWTD,  »  „ 


CQ  ”  physical  constant  evaluated  by  continuing  a  one-dimensional 

consolidation  test  past  the  end  of  primary  consolidation  for  a 
suitable  load  increment 

H  *  thickness  of  compressible  layer 

C2  ■  time  at  which  the  value  of  secondary  compression  is  desired 

tj  ■  time  at  the  beginning  of  secondary  compression;  the  time  corre¬ 
sponding  to  90  percent  of  primary  consolidation  is  sometimes  used 

For  the  silty  clay  in  this  problem,  the  secondary  settlement  is  thus  relatively 

small  compared  to  a  primary  consolidation  settlement  of  26.0  in.  If  organics 

had  been  present,  secondary  settlement  would  have  been  significantly  greater. 

Settlement  Example  2 

9.  Settlement  Example  2  illustrates  how  to  handle,  at  least  approxi¬ 
mately,  stress  distribution  in  calculating  settlement  of  stone  column  improved 
ground.  Stone  column  improved  ground  is  being  considered  as  one  design  alter¬ 
native  for  a  slightly  marginal  site  consisting  of  firm  to  stiff  sandy  silt  as 
shown  in  Figure  D4.  The  average  contact  stress  is  o  =  p/A  =  400  kips/(13  ft 

*  13  ft)  ■  2,367  psf.  The  gross  area  replacement  ratio  from  Equation  3*  is 

2 

ag  *=  Ag/A  =  (7.07  ft  )(4)/(13  ft  *  13  ft)  *  0.167  .  Now  determine  the  initial 
effective  stress  at  the  center  of  each  layer: 

Layer  1:  *  8  ft  (120  pcf)  =  960  psf 

Layer  2:  =  13  ft  (120  pcf)  +  4  ft  (123  pcf  -  62.4  pcf) 

=  1,810  psf 

Calculate  the  change  in  stress  Ao^  at  the  center  of  each  layer  using  as  an 
approximation  Boussinesq  stress  distribution  theory  for  a  square  foundation 
and  the  average  applied  stress  o  (refer  to  Lambe  and  Whitman  (1969), 

Leonards  (1962),  Sowers  (1979),  and  Winterkorn  and  Fang  (1975)  for  discussion 
of  stress  distribution  and  for  charts  and  tables  that  calculate  changes  in 
stress  due  to  foundation  loadings): 

Layer  1:  z/B  -  J-fl't  =  0,38  (02) 

A a  =  I  •  a  =  0.82  (2,367  psf)  (D3) 

z  z 

A -  1,941  psf 


D9 


Figure  D4.  Settlement  Example  2 — Rigid  Foundation 
Placed  Over  Stone  Column  Improved  Sandy  Silt 


Layer  2:  z/B  *  ^  - -  =  1.08 

La  -  I  •  a  =  0.31  (2,367  psf) 
z  z 

La  =  734  psf 
z  r 


depth  below  foundation 
foundation  width 
influence  factor 


The  change  in  stress  A calculated  above  is  the  average  stress  change  over 
the  unit  cell. 

10.  Assume  a  stress  concentration  factor  n  =  3  (an  n  value  less 
than  4  is  used  because  the  soil  is  relatively  stiff  compared  with  soft  clays) 
The  stress  ratio  in  the  sandy  silt  from  Equation  8a*  is 


1  +  (n  -  l)a 


1  +  (2.0) (0. 167) 


=  0.750 


(8a)* 


The  stress  change  in  the  sandy  silt  as  an  approximation  can  be  taken  to  equal 

U  Ao  giving  the  following  settlements  from  Equation  20*  for  layers  1  and  2: 
c  z 


5ti  “  (iro)logio (' 


960  psf  +  1941  psf  (0 
960  psf 


.750)^ 


=  1.5  in. 


(10  ft  x  12  in.)  (20)* 


»«2  - 


1810  psf  +  734  psf  (0.750)) 
1810  psf  ) 


(8  ft  x  12  in.)  (20)* 


The  total  settlement  in  the  sandy  silt  strata  is  about  =  1.9  in.  Had 
stone  columns  not  been  used,  the  settlement  would  have  been  S  =  2.4  in., 
giving  S^/S  *  1.9/2. 4  *  0.79  .  From  Equation  22*  and  Figure  C6  (19*),  St/S 
s  Uc  “  0.75  ,  which  illustrates  that  S^/S  =  is  a  quite  useful  approach 
for  preliminary  estimates  of  the  level  of  reduction  of  settlement  for  various 
stone  column  designs.  In  the  above  simplified  equation,  the  variables 
affecting  the  settlement  ratio  St/S  are  only  ag  and  n  . 

11.  Stress  distribution  can  also  be  approximately  considered  using  the 
finite  element  design  charts.  To  do  this  an  average  stress  a  is  calculated 
within  the  compressible  layer  and  used  in  the  chart  rather  than  the  stress 
actually  applied  at  the  top  of  the  layer. 

12.  In  Settlement  Example  2,  only  four  stone  columns  are  used.  Also, 
two  layers  of  sandy  silt  are  present  which  would  have  different  coefficients 

of  consolidation.  Assume  C  (and  C  )  in  one  layer  differs  from  C  (and 

v  vr  v 

C  )  in  the  other  layer  by  a  factor  of  about  2  to  5.  For  the  resulting  com¬ 
plex  three-dimensional  flow  conditions,  a  theoretically  accurate  evaluation  of 
the  time  rate  of  settlement  for  this  problem  would  be  a  major  undertaking. 

Such  a  solution  would  require  a  three-dimensional  numerical  analysis.  As  a 


rough  engineering  approximation,  however,  the  following  simplified  approach 
can  be  taken: 

a.  Consider  for  each  layer  radial  and  vertical  flow  separately  and 
use  Equation  25*  to  estimate  the  combined  results. 

b.  For  radial  flow  neglect  any  interaction  between  the  two  layers. 
Sketch  in  the  approximate  radial  flow  paths  on  a  scale  drawing 
(Figure  D5) .  Remember  that  flow  originates  from  lines  of  geo¬ 
metric  symmetry  and  moves  approximately  radially  to  the  drains. 

Consider  the  flow  to  the  drain  shown  in  the  upper  left-hand 
corner  of  Figure  D5.  An  examination  of  the  flow  paths  on  the 
figure  show  25  percent  of  the  flow  to  the  stone  column  from 
quadrant  a-o-b  is  from  infinity.  This  means  D  for  this  quad¬ 
rant  is  very  large,  and  from  Equation  28*  the  radial  time  fac¬ 
tor  =  0  .  Over  quadrants  b-o-c  and  a-o-e,  which  together 

comprises  another  25  percent  of  the  drain,  the  flow  path  length 
varies  from  infinity  at  points  b  and  a  to  short  drainage  paths 
at  points  c  and  e;  this  combined  quadrant  will  only  be  partially 


Figure  D5.  Approximate  radial  flow  paths  for  settlement 

Example  2 
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effective  in  providing  drainage.  Finally,  the  area  contributing 
flow  to  the  drain  that  lies  to  the  right  and  below  line  c-o-e 
has  short  flow  paths  that  can  be  approximated  by  an  estimated 
equivalent  unit  cell  diameter  D  =  7.5  ft  shown  in  dashed 
lines  in  the  figure.  As  an  engineering  approximation  for 
this  example,  estimate  the  time  factor  for  each  layer  using 

the  appropriate  value  of  and  =  7.5  ft.  To  crudely 

consider  the  effectively  is  very  large  over  between  25  to 

50  percent  of  the  drain,  reduce  the  time  factor  by  about 
(25  percent  +  50  percent)/2  =  37.5  percent  or  0.4  (multiply  the 
calculated  time  factor  T  by  0.6  or  use  0.5  to  be  a  little 
more  conservative). 

To  illustrate  this  approximate  approach  assume  for  Layer  1 

2 

C,_  =  5C  and  C  =  0.2  ft  /day.  Then  at  the  end  of  2  months 
vr  v  v 

the  radial  time  factor  would  be  estimated  from  Equation  28*  as 

Tr  "  Cvr  *  t/De  =  5(0,2  ft2/dayX2  *  31  days)/(7.5  ft)2  =  1.10  . 

Reducing  the  time  factor  to  approximately  consider  partial 
drainage  gives  T^  =  0.5  (1.10)  =0.55  .  Assume  the  stone 

column  diameter  is  effectively  reduced  by  one-fifth  to  account 
for  smear,  giving,  from  Equation  29a*,  n*eqUiv  “  2,3  ft/ (3  ft 

x  0.2)  =  12.5.  Then  from  Figure  D3  (43*)  the  degree  of  radial 
consolidation  =  0.91  .  Conservatively  neglecting  vertical 

drainage  in  Layer  1,  the  settlement  after  2  months  of  Layer  1 
is  Stl  s  1.5  in.  (0.9)  =  1.35  in.  As  would  be  expected, 

consolidation  occurs  rapidly  in  the  sandy  silt. 


* 

r 

w 


r 

e 


* 

nequiv 


(29)* 


where 


r* 

w 

r 

e 

n* 

equiv 


radius  of  drain  without  smear 
radius  of  unit  cell 

radius  of  unit  cell  divided  by  equivalent  radius 
of  drain  without  smear 


Since  in  this  example  =  5  for  Layer  1,  vertical  com¬ 

pared  to  radial  consolidation  would  be  relatively  slow  and  was 
conservatively  neglected.  However,  if  the  effect  of  vertical 
drainage  on  the  time  rate  of  consolidation  is  desired,  the 
presence  of  two  layers  greatly  complicates  vertical  time  rate 
of  consolidation  computations.  If  of  the  more  permeable 

layer  is  more  than  20  times  of  the  less  permeable  layer 

the  following  simplified  approach  can  be  used  (Larabe  and 
Whitman  1969):  (1)  assume  consolidation  occurs  in  two  stages; 


D13 


(2)  in  Stage  1,  calculate  consolidation  f  the  more  permeable 
layer,  assuming  no  drainage  at  the  interface  between  the  two; 
and  (3)  in  Stage  2,  calculate  consolidation  in  the  least  per¬ 
meable  layer,  assuming  drainage  at  the  interface.  If  C  of 
one  layer  is  less  than  20  times  C  of  the  other,  the  v 
approximate  method  described  in  NA^FAC  DM-7  (Department  of  the 
Navy  1971)  can  be  followed  or  numerical  methods  can  be  used 
(Lambe  and  Whitman  1969) . 

The  above  methods  are,  of  course,  quite  crude  and  should  be 
considered  "ball  park"  in  accuracy.  They  do  give  a  rational 
way  of  approaching  a  very  complicated,  three-dimensional  time 
rate  of  primary  consolidation  problem. 
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APPENDIX  E:  EXAMPLE  STABILITY  PROBLEM 


1.  This  example  illustrates  how  to  handle  the  geometric  and  material 
parameters  required  for  setting  up  a  slope  stability  problem  for  analysis 
using  the  the  Profile  Method  described  in  Chapter  III,  Barksdale  and  Bachus 
(I983t). 

2.  A  I5-ft  high  embankment  is  to  be  placed  over  a  soft  clay  as  illus¬ 
trated  in  Figure  El.  Because  of  the  low  shear  strength  of  the  soft  clay,  use 
a  stress  concentration  factor  n  of  2.0,  and  an  angle  of  internal  friction 

$  of  the  stone  column  of  42°.  The  saturated  unit  weight  of  the  stone  is 
s 

125  pcf.  For  the  first  trial  design,  use  five  rows  of  stone  columns  laid  out 
as  shown  in  Figure  El.  An  equilateral  triangular  grid  will  be  used  having  a 
trial  spacing  s  *  6.5  ft.  The  stone  column  diameter  is  estimated  to  be 
3.5  ft,  giving  an  area  replacement  ratio  of 


a 

s 


0.263 


(5b)* 


3.  The  plan  view  of  the  stone  column  grid  used  to  improve  the  site  is 
shown  in  Figure  El(b).  As  shown  in  the  figure,  stone  columns  replace  only 
26  percent  of  the  total  volume  of  the  soft  clay  (i.e.,  ag  =  0.263).  Further, 
in  performing  a  conventional  stability  analysis,  the  materials  are  assumed  to 
extend  for  an  infinite  distance  in  the  direction  of  the  embankment.  Typi¬ 
cally,  the  analysis  is  then  performed  on  a  1-ft-wide  slice  of  embankment.  To 
use  the  profile  method  the  discrete  stone  columns  must  therefore  be  converted 
into  equivalent  stone  column  strips  extending  along  the  full  length  of  the 
embankment  as  follows: 


A 

s 


3.14 


(3.5  ft)2 
4 


9.62  ft 


(El) 


Note:  Equation,  table,  or  figure  numbers  with  an  asterisk  refer  to  the  Fed¬ 

eral  Highway  Administration  stone  column  report  (Barksdale  and  Bachus, 
1983). 

t  References  cited  in  this  Appendix  can  be  found  on  the  last  page. 
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Figure  El.  Stability  Example  Problem  1 

The  length  tributary  to  each  stone  column  in  the  direction  of  the  embankment 
is  s  m  6.5  ft.  Therefore,  a  solid  strip  having  the  same  area  and  volume  of 
stone  would  have  a  width  w  of 

A  ? 

_s  _  9.62  ft 

s  6.5ft 


w 


1.48  ft 


(E 


2 

Tha  total  width  of  the  tributary  area  equals  A  /(a  s)  =  9.62  ft  /  (0.263 

s  s 

x  6.5  ft)  *  5.63  ft,  which  is  the  stone  column  spacing  0.866s  in  the 
direction  perpendicular  to  the  embankment  length. 

4.  Now  determine  the  characteristics  of  fictitious  strips  (Figure  E2) 
that  must  be  added  to  handle  the  effect  of  stone  column  stress  concentration 
in  the  stability  analysis.  Let  the  thickness  T  of  a  fictitious  strip  be 
0.3  ft  under  the  full  embankment.  Note  that  in  this  example  no  stone  columns 
are  actually  used  under  the  full  embankment  height  for  the  first  trial.  How¬ 
ever,  stone  column  row  5  is  located  so  that  the  edge  of  the  tributary  area  is 
just  at  the  break  in  the  embankment.  The  unit  weights  calculated  for  the  full 
embankment  height  can  be  used  for  each  strip,  with  the  thickness  of  the  strips 
varying  from  zero  at  the  toe  to  0.3  ft  at  the  break  (Figure  E2) .  An  examina¬ 
tion  of  Equations  33*  and  34*  shows  that  this  method  gives  the  proper  stress 
concentration  in  each  strip.  The  thickness  of  the  boundaries  T  of  each 
zone  is  calculated  in  Table  El. 


r 


Figure  E2.  Zones  used  for  computer  idealization  for  Stability  Example 

Problem  1 


Note:  -  Thickness  of  fictitious  layer  at  location  shown  in  Figure  E2.  z 


is  defined  as  the  horizontal  distance  from  the  break  in  the  embank¬ 
ment  slope  to  the  point  at  which  the  thickness  T  is  required.  Then 


from  similar  triangles  T^/OO  ft  -  z)  *  0.3  ft/30  ft. 


5.  The  unit  weights  to  use  in  the  fictitious  strips  are  calculated  as 
follows : 

n  2.0  , 


s  [1  +  (n  -  1)  a  ]  [1  +  (0.263)] 


_ 1 _  _ L0 _ n  7Q? 

V  [1  +  (n  -  1)  a  ]  “  [1  +  (0.263)]  "  *  y 

s 


6.  The  correct  unit  weight  to  use  above  each  stone  column  in  the 
fictitious  strip  is 


Yj  -  (,,  -  DYjHVT  -  ‘).0iP.pef)»5  ft) 

■  3,480  pcf 


s 


fiti 

&  I 


KM,« 

iViy 


and  the  unit  weight  to  use  above  the  soil  in  each  fictitious  strip  is 


where 


Yf  -  (uc  -  DYjH'/T 
*  -1,248  pcf 


(0.792  -  1) (120  pcf ) (15  ft) 
0.3  ft 


(34)* 


*  unit  weight  of  embankment 
H'  *  height  of  embankment 

7.  Material  properties  and  zones  are  as  follows  (refer  to  Figure  E2) : 


Zone  1:  yw  “  120  pcf,  c  =  50  psf,  <p  =  28° 

Zone  2:  y  =  0,  c  *=  0  ,  <f>  -  0 
w 

Zones  3,  5,  7,  9,  11,  13:  y  =  -1,248  pcf,  $  =  0  ,  c  =  0 

Zones  4,  6,  8,  10,  12:  y  =  +3,480  pcf,  <p  =  0  ,  c  =  0 

Zones  14,  16,  18,  20,  22,  24:  y  =  110  pcf,  <f>  =  0  ,  c  =  350  psf 

Sa  L 

Zones  15,  17,  19,  21,  23:  y  =  125  pcf,  $  =  42°  ,  c  =  0 

Sol. 

8.  The  calculated  safety  factor  of  the  slope  is  shown  in  the  table 
below  for  the  following  conditions:  (a)  no  improvement  with  stone  columns, 

(b)  the  stone  column  improvement  shown  in  Figure  El  using  a  stress  concentra¬ 
tion  factor  n  *  1  ,  and  (c)  the  same  level  of  improvement  with  n  =  2.0  (the 
critical  circle  for  this  condition  is  shown  in  Figure  El).  A  simplified 
Bishop  analysis  was  performed  using  the  GTICES  Lease  II  computer  program  pre¬ 
pared  by  A.  W.  Dawson  (1972). 


Table  E2 
Safety  Factors 


Coordinate  ^ 


Case 

n 

X 

ft 

y 

ft 

R 

ft 

S.F. 

Comment 

1. 

No  S.C.(2;> 

_ 

14.20 

27.00 

43.00 

1.07 

Base  Failure 

2. 

S.C. 

1 

2.90 

26.00 

42.00 

1  .38 

Base  Failure 

3. 

S.C. 

2.0 

2.90 

26.00 

34.75 

1.65 

See  Figure  El 

Notes:  1.  Coordinates  of  critical  circle  (refer  to  Figure  El  for  location  of 

x  and  y  axles) . 

2.  Notation:  S.C.  =  stone  column;  S.F.  =  Safety  factor;  R  =  radius 
of  critical  circle. 
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APPENDIX  F:  EXAMPLE  PROBLEM  OF  PORE  PRESSURE  DISSIPATION 
DURING  AN  EARTHQUAKE 


1.  Stone  columns  approximately  30  ft  in  length  are  to  be  used  in  loose 
sands  to  (a)  densify  the  sand,  (b)  provide  additional  lateral  shear  strength, 
and  (c)  dissipate  pore  pressures  generated  during  an  earthquake.  The  site 
conditions  consist  of  alluvial,  generally  loose,  fine-to-medium-grained  sands 
interbedded  with  silt  and  clay  layers  and  lenses.  The  maximum  triangular  pat¬ 
tern  stone  column  spacing  s  is  required  to  prevent  liquefaction  from  occur¬ 
ring  within  the  sand  layers  during  a  design  earthquake. 

2.  An  equilateral  triangular  pattern  of  stone  columns  is  to  be  employed. 

The  stone  columns  will  be  constructed  using  a  processed  river  gravel  having  a 

maximum  size  of  1.5  in.,  and  a  minimum  size  of  0.375  to  0.5  in.  The  completed 

diameter  of  the  stone  columns  is  estimated  to  be  3.0  ft.  The  sands  have  a 

horizontal  permeability  of  about  0.007  ft/min,  and  a  coefficient  of  compress- 

_  6  2 

ibility  m^  of  2.5  *  10  ft  /lb.  Both  permeability  and  compressibility 
were  estimated  for  conditions  after  stone  column  construction.  The  character¬ 
istics  of  the  earthquake  as  determined  from  a  separate  detailed  earthquake 
analysis  are  as  follows: 


N 

N. 


1.5 


t,  =  0.5  min 
a 


(FI) 

(F2) 


where 

N  =  number  of  uniform  stress  cycles  equivalent  to  earthquake 
eq 

N^  *  number  of  cycles  to  liquefaction 

t,  =  time  over  which  the  N  uniform  stress  cycles  of  the  earthquake 
d  14  j  eq 

are  applied 

3.  The  stone  columns  are  to  be  constructed  using  the  pull-down  pro¬ 
cess.  The  jets  on  the  outside  of  the  vibroflot  (casing)  and  the  horizontal 
vibrator  will  both  be  operated  during  construction.  Use  of  this  construction 
technique  will  prevent  sand  from  being  deposited  within  the  open-graded  stone 
backfill.  Although  the  effects  of  smear  should  be  relatively  small  using  this 
technique,  some  smear  in  the  sand  may  result  due  to  dragging  of  soil  from  the 
cohesive  layers.  Therefore,  assume  the  effective  stone  column  diameter  is 
one-half  of  the  actual  diameter  of  3  ft. 

FI 
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4.  The  permeability  of  the  stone  column  backfill  material  as  measured 
in  the  laboratory  is  35  ft/min.  Hence  the  ratio  of  the  vertical  stone  column 
drain  to  the  in  situ  sand  permeability  is  35  ft/min/0.007  ft/min  =  5,000  . 

This  value  is  considerably  greater  than  the  ratio  of  200  indicated  by  Seed  and 
Booker  (1976) t  for  water  to  readily  move  vertically  through  the  stone  column 
during  an  earthquake. 

5.  To  determine  the  required  stone  column  spacing,  first  calculate  the 
dimensionless  time  factor 


4k  t , 
r  d 


4(0.007  ft/min)  (0.50  min) 

1 /^3.  /o  c  ..  ,  n"”^  c  I  ■ ,  c  \ 


Y  m  V  (62.4  lb/ft3)  (2.5  x  10~6  ft2/lb)  (1.5  ft)2 
w  v3 


T  .  =  39.9 
ad 


where 


k^  =  horizontal  permeability  of  the  sand  stratum 

Yw  =  unit  weight  of  water 
m^  =  coefficient  of  compressibility 
D  =  diameter  of  the  stone  column 

Note  that  the  actual  stone  column  diameter  of  3  ft  has  been  reduced  by 
one-half  in  the  above  equation  to  crudely  account  for  possible  smear. 

6.  To  be  conservative  use  a  safety  factor  of  1.5  with  respect  to  the 
greatest  pore  pressure  ratio  in  the  sand,  i.e.,  let  the  greatest  allowable 
value  of  the  pore  pressure  ratio  r^  =  1/SF  =  1.0/1. 5  =  0.67  .  Since  Neq/N£ 
=  1.5,  interpolation  using  Figures  3,  4,  and  5  (see  main  text)  is  required  as 


follows : 


Figure 


N  /N 


Required 

D/D 


By  interpolating  from  a  plot  of  the  above  data  for  ^^/N^  =  1.5  ,  the 
required  ratio  D/D  =  0.16  .  Hence  Dg  =  D/0.16  =  1.5  ft/0.16  =  9.4  .  Since 
the  stone  column  spacing  s  =  D  /1.05  ,  then  s  =  9.4  ft/1.05  =  8.95  ft  .  A 

t  References  cited  in  this  Appendix  can  be  found  on  the  last  page. 
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9-ft  triangular  spacing  should  therefore  provide  effective  lateral  dissipation 
of  excess  pore  pressure  during  the  design  earthquake.  Of  course,  all  other 
pertinent  factors  such  as  required  sand  relative  density,  settlement,  and 
shear  should  be  considered  in  selecting  a  final  stone  column  spacing. 
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APPENDIX  G:  NOTATION* 


a 

s 

A 

A 

s 

B 

c 

C 

c 


C 


a 


C 

v 


C 

vr 


D 


D 

e 


e 

o 

E 


E 

c 

E 

8 

F 


h 

H 


H' 


Area  replacement  ratio,  Ag/A 
Total  area  within  the  unit  cell 
Stone  column  area 
Foundation  width 
Cohesion  of  soil 

Virgin  compression  index  of  cohesive  soil  from  one-dimensional  con¬ 
solidation  test 

Coefficient  of  secondary  compression,  C^  *  AH/  [■  log10 
Coefficient  of  consolidation  in  vertical  direction  (Equation  27) 
Coefficient  of  consolidation  in  radial  direction  (Equation  28J 
Constructed  diameter  of  stone  column  (Figures  13  and  14) 

Equivalent  diameter  of  unit  cell  (Equations  1  and  2) 

Initial  void  ratio  of  cohesive  soil 
Modulus  of  elasticity 

Modulus  of  elasticity  of  thin  boundary  around  the  unit  cell  used 
in  nonlinear  finite  element  analysis 

Modulus  of  elasticity  of  soil  within  the  unit  cell 

Modulus  of  elasticity  of  the  stone  column 

Shear  force  on  upper  failure  surface  in  stone  column  undergoing 
local  bearing  failure  (Appendix  B) 

Vesic  cavity  expansion  factors  (Figure  16) 

Depth  of  fill  above  the  stone  column 

Vertical  height  (or  increment)  of  stone  column  treated  ground  over 
which  settlements  are  calculated 

Height  of  embankment  in  stability  analysis  (Figure  46) 


Note:  Equation,  table,  or  figure  numbers  refer  to  the  Federal  Highway 
Administration  stone  column  report  (Barksdale  and  Bachus,  1983). 

*  References  cited  in  this  Appendix  can  be  found  in  the  References  at  the 
end  of  the  main  text. 
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Rigidity  index  used  in  Vesic  cavity  expansion  theory  (Equation  13) 
Horizontal  permeability  of  the  sand  stratum 
Permeability  of  soil  in  radial  direction  (Figure  45) 

Permeability  of  smear  zone  in  radial  direction  (Figure  45) 

Permeability  of  soil  in  vertical  direction 

Coefficient  of  at-rest  earth  pressure 

Length  of  stone  column 

Coefficient  of  compressibility 

Stress  concentration  factor,  a  /a  (Figure  14) 

s  c 

Ratio  of  the  unit  cell  radius  to  the  radius  of  the  drain  (stone 
column  radius  less  smear  zone  thickness),  n*  =  re/rw  (Figure  43) 

Equivalent  value  of  n*  for  a  drain  without  smear,  n*  =  r  /r* 
(Figure  44)  6q  6  W 

Number  of  drainage  surfaces  at  the  top  and  bottom  of  the  layer 
(N  =  1  or  2) ;  also  normal  force  on  lower  failure  surface  in  stone 
column  undergoing  local  bearing  failure  (Appendix  B) 

Ultimate  bearing  capacity  factor  of  stone  column  (Equation  50) 

Number  of  uniform  stress  cycles  equivalent  to  earthquake 

Number  of  cycles  to  liquefaction 

Ultimate  lateral  resistance  of  clay  acting  on  critical  wedge  for  a 
local  bearing  failure  of  stone  column  (Appendix  B) 

Ultimate  loading  that  can  be  applied  over  the  unit  cell  area 

Mean  isotropic  stress,  q  =  ^  +  o^)/3 

Ultimate  bearing  capacity 

Ultimate  bearing  capacity  of  stone  column 

Radius  of  the  unit  cell  (Figure  45) 

Greatest  pore  pressure  ratio 

Greatest  average  pore  pressure  ratio 

Radius  of  smear  zone  (Figure  45) 
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r  Radius  of  the  drain  usually  taken  as  the  radius  of  the  stone  column 
less  the  thickness  of  the  smear  zone  (Figure  45) 

r*  Radius  of  equivalent  drain  without  smear  (Figure  45) 

R  Radius  of  the  stone  column 

s  Center  to  center  spacing  of  stone  columns  (Figure  13) 

s*  Ratio  of  the  radius  of  smear  zone  to  radius  of  the  drain, 
s*  =  r  /r  (Figure  44) 

8  W 

S  Settlement  of  unimproved  ground 

S  Settlement  occurring  in  an  increment  H  of  stone  column  treated 
ground 

S^  Primary  consolidation  settlement  at  time  t 

Settlement  of  a  single  stone  column  (Figure  50) 
t  Time 


t.  Time  over  which  the  N  uniform  stress  cycles  of  the  earthquake  are 
d  ,  eq 

applied 


T  Shear  force  on  lower  failure  surface  in  stone  column  undergoing 
local  bearing  failure  (Appendix  B) 


T  Assumed  thickness  of  fictitious  strip  of  soil  used  to  obtain  proper 
stress  concentration  in  a  computer  stability  analysis  (Figure  46) 

2 

T  Time  factor  for  radial  drainage,  T  =  C  t/(D  )  (Figure  43) 
r  r  vr  e 

2 

T  Time  factor  for  vertical  drainage,  T  =  C  t/ (H/N)  (Figure  42) 
z  Z  V 

T  ,  Dimensionless  time  factor 
ad 

u  Induced  pore  pressure 

U  Average  degree  of  consolidation  considering  both  vertical  and 
radial  drainage,  U  =  1  -  (1  -  U  )(1  ~ 

U  Average  degree  of  consolidation  in  radial  (horizontal)  direction 
r  (Figure  43) 

U  Average  degree  of  consolidation  in  vertical  direction  (Figure  42) 
z 

w  Width  of  equivalent,  continuous  stone  strip  used  in  a  stability 
analysis  w  *=  A^/s  (Figure  46) 

W  Effective  normal  force  exerted  on  upper  failure  wedge-local  bearing 
failure  (Appendix  B) 
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Effective  weight  of  stone  in  failure  wedge-local  bearing  failure  of 
stone  column  (Appendix  B) 

Effective  vertical  force  exerted  on  the  circular  arc  failure  sur¬ 
face  or  the  upper  surface  of  the  failure  wedge  for  local  bearing 
failure  (Appendix  B) 

Depth  below  ground  surface 

Depth  of  circular  arc  failure  below  top  of  stone 

Inclination  of  lower  failure  surface  in  a  stone  column  undergoing  a 
local  bearing  failure  (Appendix  B) 

A  parameter  characterizing  the  shape  of  the  pore  pressure  genera¬ 
tion  curve 

Inclination  of  shear  surface  with  respect  to  the  horizontal 

Saturated  (wet)  unit  weight  of  cohesive  soil 

Weight  of  fictitious  soil  strip  for  use  in  computer  stability 
analysis,  y^  =  (pc  -  Oy^H/T  (Figure  46) 

Weight  of  fictitious  soil  strip  for  use  in  computer  stability 
analysis,  y®  =  (pg  -  Oy^H/T  (Figure  46) 

Saturated  (wet)  unit  weight  of  stone  column 

Bouyant  unit  weight  of  stone  in  failure  wedge-local  bearing  failure 
Unit  weight  of  water 

Average  unit  weight  of  material  within  unit  cell 

Unit  weight  of  embankment  in  stability  analysis  (Figure  46) 

Reduction  factor  for  local  bearing  failure  of  a  stone  column 
(Appendix  B) 

Reduction  factor  to  apply  to  measured  field  vane  shear  strengths 
(Figure  73) 

Ratio  of  st' cas  Ln  cohesive  soil  to  average  stress,  y  =  a 
*  o^/a  (Equation  8a)  ~  C 

Ratio  of  stress  in  stone  column  to  average  stress,  y  =  a  /o 
(Equation  8b) 

Poisson’s  ratio 


Poisson's  ratio  of  soil 
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Poisson's  ratio  of  stone  column 

Average  stress  acting  over  the  unit  cell  area  due  to  the  applied 

loading  (Figure  14) ;  stress  distribution  should  be  considered  with  | 

depth  where  important  (Figure  40) 

Average  stress  acting  over  the  soil  in  the  unit  cell  (Figure  14)  j 

Average  stress  acting  over  the  stone  column  (Figure  14)  j 

Initial  effective  stress  in  cohesive  soil  before  stone  column 

construction  j 

Average  of  initial  and  final  stress  state  applied  to  the  cohesive  ! 

soil;  o  is  used  in  Equation  47  to  calculate  E  for  consolidation  1 

test  results  ! 

Major  principal  stress 
Minor  principal  stress 

Shear  strength  in  cohesive  soil  on  failure  surface  in  a  stability 
analysis 

Shear  strength  in  stone  on  failure  surface  in  a  stability  analysis 

Earthquake  equivalent  uniform  stress 

Angle  of  internal  friction  of  soil 

Angle  of  internal  friction  of  stone  column 

tan  o  -  tan  6  (Equation  56) 


Figure  Gl.  Unit  cell  idealization 


